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^|- . ABSTRACT 

in 



Far-infrared continuum data from the COBEj DIRBE instrument were com- 



O 

bined with Nagoya 4-m 13 C0 J = 1 — > spectral line data to infer the 
1 multiparsec-scale physical conditions in the Orion A and B molecular clouds, 

using 140 /zm/240 /im dust color temperatures and the 240/xm/ 13 CO J = 1 — > 
intensity ratios. In theory, the ratio of far-IR, submillimeter, or millimeter con- 
tinuum to that of a 13 CO (or C 18 0) rotational line can place reliable upper limits 
on the temperature of the dust and molecular gas on multi-parsec scales; on 
such scales, both the line and continuum emission are optically thin, resulting 
in a continuum-to-line ratio that suffers no loss of temperature sensitivity in the 
high-temperature limit as occurs for ratios of CO rotational lines or ratios of 
continuum emission in different wavelength bands. 

Two-component models fit the Orion data best, where one has a fixed- 
temperature and the other has a spatially varying temperature. The former 
represents gas and dust towards the surface of the clouds that are heated pri- 
marily by a very large-scale (i.e. ~ 1 kpc) interstellar radiation field. The latter 
represents gas and dust at greater depths into the clouds and are shielded from 
this interstellar radiation field and heated by local stars. The inferred phys- 
ical conditions are consistent with those determined from previously observed 
maps of 12 CO J = 1 — > and J = 2 — > 1 that cover the entire Orion A and B 
molecular clouds. The models require that the dust-gas temperature difference 
is 0±2 K. If this surprising result applies to much of the Galactic ISM, except in 
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unusual regions such as the Galactic Center, then there are a number implica- 
tions. These include dust-gas thermal coupling that is commonly factors of 5 to 
10 stronger than previously believed, Galactic-scale molecular gas temperatures 
closer to 20 K than to 10 K, an improved explanation for the N(H 2 ) /I(CO) con- 
version factor (a full discussion of this is deferred to a later paper), and ruling 
out at least one dust grain alignment mechanism. The simplest interpretation 
of the models suggests that about 40-50% of the Orion clouds are in the form 
of cold (i.e. ~ 3-10 K) dust and gas, although alternative explanations are not 
ruled out. These alternatives include the contribution to the 240 /mi continuum 
by dust associated with atomic hydrogen and reduced 13 CO abundance towards 
the clouds' edges. Even considering these alternatives, it is still likely that cold 
material with temperatures of ~ 7-10 K still exists. If this cold gas and dust are 
common in the Galaxy, then mass estimates of the Galactic ISM must be revised 
upwards by up to 60%. 

The feasibility of submillimeter or millimeter continuum to 13 CO line ratios 
constraining estimates of dust and molecular gas temperatures was tested. The 
model fits allowed the simulation of the necessary millimeter-continuum and 
13 CO J = 1 — ► maps used in the test. In certain "hot spots" — that have 
continuum-to-line ratios above some threshold value — the millimeter contin- 
uum to 13 CO ratio can estimate the dust temperature to within a factor of 2 
over large ranges of physical conditions. Nevertheless, supplemental observations 
of the 13 CO J = 2 — > 1 line or of shorter wavelength continuum are advisable 
in placing lower limits on the estimated temperature. Even without such sup- 
plemental observations, this test shows that the continuum-to-line ratio places 
reliable upper limits on the temperature. 

Subject headings: ISM: molecules and dust — Orion 



Introduction 



While interesting in themselves, molecular clouds provide insights into star forma- 
tion. Since stars form in and from molecular clouds, knowing the physical conditions 
within these cloud s is essentia l for a complete understanding of star formation. As men- 
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the interstellar medium in general, can also be observed in the millimeter, submillimeter, and 
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19881 ). Estimating physical parameters like tem- 



perature, and sometimes density, requires using the ratios of intensities of spectral lines or 
of the continuum at different wavelengths. Given that each of these ratios is dependent on 
the ratio of two Planck functions at two different wavelengths, they often lose temperature 
sensitivity at higher temperatures. While there are methods of addressing this shortcoming, 
having two tracers of molecular gas with different dependences on the temperature would 
complement other methods of tracing warm dust or molecular gas. This is especially true if 
the tracers are optically thin, because low opacity emission is more sensitive to the physical 
parameters of the bulk of the gas, rather than in just the surface layers. 

One such pair of tracers is a rotational line of an isotopologue of CO, such as that of 13 CO 
or C 18 0, and the submillimeter continuum. Both of these tracers are optica lly thin on the 
scales of m any parsecs, w hich are the scales of interest for the current work. ISchloerb et al. 



(119871 ) and ISwartz et al.l (Il989l ) showed that the intensity ratio of an optically thin isotopic 
CO line emission to submillimete r continuum emission can estimate the temperature of gas 
and dust in molecular clouds. The lSchloerb et al.l (119871 ) expression for this ratio goes roughly 
like T 2 in the high-temperature limit. Accordingly, the I I/ (submm)/I(C 18 J = 1 — ► 0) 
and I^(submm)/I( 13 CO J = 1 — > 0) ratios are actually more sensitive to temperature as 
that temperature increases. This is in stark contrast to ratios of rotational lines of a given 
isotopologue of CO and to ratios of continuum intensities at different frequencies, which 
lose sensitivity to temperature in the high-temperature limit. The I^(submm)/I( 13 CO) ratio 
can then serve as the needed diagnostic of high gas/dust temperatures, provided that the 
shortcomings and complications of these tracers can be overcome or at least mitigated. These 
complications include variations in the 13 CO-to-dust mass ratio, non-molecular phases of the 
ISM along the line of sight, variations of gas density, variations in dust grain properties, 
appreciable optical depth variations in the 13 CO line used, and others (see the Introduction 
of Paper I for more details). Such complications are often reduced in the case of observations 
on multi-parsec scales, because spatial gradients on such scales are generally smaller than 
the extremes that occur on very small scales. Consequently, testing the reliability of the 
Iy(submm)/I( 13 CO) ratio as high-temperature diagnostic is best carried out with observations 
of a molecular cloud, or of clouds, on multi-parsec scales. 
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The Orion A and B molecular clouds were chosen as the clouds for testing the 
I ;/ (submm)/I( 13 CO) ratio 's diagnostic abilit y. They have been mapped in the 13 CO J 



1 — ► line (see, e.g.. iNagahama et al 



19981 ) and in the far-IR by IRAS flBallv et al.lll991l ) 



and COBE/DIRBE f lWall et al.l Il996l . W96 hereafter). Avoiding the complication of the 
emission of stochasticall y heated dust grai ns requires far-IR observations at wavelengths 
longer than 100 jum (e.g. iDesert et al.l ll990L W96). Accordingly, the far-IR observations of 
COBE/DIRBE were used instead of IRAS because the former has two bands — A = 140 fim 
and 240 /im - - longward of 100 /im, whereas the latter does not. The Orion A and B 
clouds were chosen for this study because they have the advantages that they are bright in 
13 CO J = 1 —>■ and at far-IR wavelengths, are out of the Galactic plane to avoid confusion 
with foreground and background emission, are several degrees in size so as to accommodate 
many DIRBE beams, and have the best ran ge of dust temperatures at the DIRBE resolution 



of 0?7 (see the Introduction of Paper I and \COBEl DIRBE Explanatory Supplement! 11998 



for more details). Therefore, the 1^(240 /xm)/I( CO) ratio, hereafter called r 240 , was plotted 
against the 140 /zm/240 /im dust color temperature, or Td c , to test the r 240 's ability to recover 
molecular cloud physical conditions. Physical models were applied to these data and physical 
conditions were inferred in Paper I. Th e reliability of the model results were tested with 
simulated data in Paper II (]Wallll2007al ). 

The next section summarizes the model results (i.e. Paper I) and the results of the 
simulations (i.e. Paper II). Section [3] then discusses general systematic effects that had not 
been treated previously. Section @] gives the scientific implications of the results. 



2. Review of the Results of the Modelling and of the Simulations 

The details of the treatment of the data and of the modeling and its results are found in 
Paper I. After a subtraction of large-scale emission from the Orion 140 /im and 240 /im maps 
representing foreground/background emission not associated with the Orion clouds (such 
subtraction was not necessary and, therefore, not applied to the 13 CO J = 1 — ► map), one- 
component and two-component model curves were fitted to the observational data in the 
r 240 versus Td c plot. There were two types of one-component models: LTE and LVG (a type 
of non-LTE model). There were also two types of two-component models (both types being 
using the LVG code): simple two-component models and two-subsample, two-component 
models. These models all adopted some form of the following assumption: 

The only physical parameters that change from one line of sight to the next are the dust tem- 
perature, T d , and the gas kinetic temperature, T K , while maintaining a constant difference, 
AT = T d — T K . Other physical parameters such as gas density, dust-to-gas mass ratio, dust 
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mass absorption coefficient, etcetera are assumed to be constant from position to position. 

This is referred to as the basic assumption. In the case of the one-component, LTE models, 
this means that the only fitted parameter was AT, while the Td and T K freely varied from 
position to position. For the one-component, LVG models, the fitted parameters were AT, 
the 13 CO column density per velocity interval, N( 13 CO)/Av, and the molecular hydrogen 
density, n(H 2 ), while the T d and T K freely varied from position to position. Put very explic- 
itly, the basic assumption applied to the one-component, LVG models means that the AT, 
N( 13 CO)/Av, and n(H2) were assumed to be spatially unchanging and therefore are the pa- 
rameters to be determined from the model fits. For the simple two-component models, there 
was a component 0, representing dust and gas in the surface of the clouds and largely heated 
by a large-scale interstellar radiation field (ISRF), and a component 1, representing dust and 
gas deeper into the clouds heated by local stars and a large-scale ISRF attenuated by the 
surface layers of gas and dust. The physical parameters of component were spatially un- 
changing and the physical parameters of component 1 were also spatially unchanging, except 
for T d and T K . The component-0 parameters were the dust temperature, T d0 , the column 
density per velocity interval, Nc °^ v C °' > ; the density, n c o, and the filling factor relative to com- 
ponent 1, Co. The component-1 parameters were Ncl ^ v C °^ and n c \. (The component-1 dust 
temperature, T dl , varied so as to generate a curve in the r 240 versus T dc plot. Component 
0, for example, would only generate a single point in this plot if there were no component 1 
contributing to the model output.) One more parameter derived from the model fit was 
the AT, which was assumed to be the same for both components. The two-subsample, two- 
component models were similar to the simple two-component models, except that the models 
were fitted to two separate subsamples within the sample of data points: the points with 
T dc < 20 K and those with T dc > 20 K. Having two subsamples allowed a better fit to the 
T dc > 20 K points; the fits are normally dominated by the T dc < 20 K subsample of points, 
often preventing good fits to the T dc > 20 K points. The resultant parameter values for all 
the model fits are summarized in Table [TJ (Notice that the two-component model results for 
the T dc > 20 K subsample are shown for the two-subsample, two-component models, whereas 
in Table 2 of Paper I the one-component model results were shown for this subsample.) 

To check the results, the systematics were tested. This was done by applying scale 
factors to the model curves that represented the effect of systematic uncertainties. These 
were the uncertainties most directly related to the comparison between the model curves 
and the observational data: the calibration of the observed 1^(240 /im)/I( 13 CO) ratio, the 
uncertainty in the dipole moment of CO, the uncertainty in the 13 CO abundance, and the 
uncertainty in the dust optical depth to total gas column density. A very rough uncertainty of 
20% was adopted for each of these uncertainties (see Paper I for details). These uncertainties 
are independent and, when added in quadrature, give a total systematic uncertainty of 40% 
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for the ratio of the model curve to the observed data. Therefore, the scale factors applied to 
the model curves ranged between 0.6 and 1.4. They were chosen to change in steps of 0.2; 
the scale factors used were 0.6, 0.8, 1.0, 1.2, and 1.4. Each scale factor was multiplied by the 
model curve before each fit for all the model types. This gave a range of fitted values for each 
parameter and this range represents the systematic uncertainty for the given parameter. 

For the two-component models an additional test was to slightly shift the starting search 
grid before running each fit. This also gave a range of results that was comparable to the 
range found from changing the scale factor. 

As a further check on the results, masses and beam filling factors were derived for each 
of the model types. Specifically, the gas-derived column densities were compared against 
the dust-derived column densities as a self-consistency check: if the model curve acceptably 
fit the data in the r 240 versus T dc plot, then there should be little scatter in the column 
density versus column density plot. As an additional check, the beam filling factors should 
be physically meaningful; they should be < 1, given that they are area filling factors. 

One flaw of the above-mentioned tests is that they cannot guarantee that the true 
values are within the ranges of results found from changing the scale factor or the starting 
grid. The method employed here could be biased to ranges of results that are far from the 
correct values. Consequently, a third series of tests was performed by fitting the one- and 
two-component models to simulated data (see Paper II). The results of the modeling and 
the tests are listed in Table 2 of Paper II. (See Table [2] of the current paper for an updated 
version of that table.) 

The most basic result of this modeling was that the two-component models fit the data 
better than the one-component models at the 99.9% confidence level, according to the F- 
test. The dust-gas temperature difference, AT, was found to be zero to within 1 or 2 K. 
The component-0 dust temperature, Tdo, was found to be 18 K (uncertainty to be discussed 
in Section [3]). The other parameters, such as the column densities per velocity interval and 
volume densities for the two components, were much less certain. This is understandable 
given that the 13 CO J = 1 —>■ line emission is well approximated by the optically thin, LTE 
limit for much of the gas of the Orion clouds. Consequently, only rough lower limits could be 
applied to the densities and rough upper limits to the column densities per velocity interval. 
The lower limit on the column density per velocity interval of component 1 is simply the 
column density per velocity interval of the clouds on the scale of the DIRBE beam. This 
lower limit for component is near the lower limit of the master search grid. 

One important consequence of the two-component models is that there is about 60% 
more mass than would be inferred from the simpler one-component models. This extra 
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material is in the form of cold (i.e. ~ 3-10 K) dust and gas. 

3. Considerations of Systematics 

In this section we consider the systematics that are less directly connected with the 
comparison between models and data. Specifically, we examine the effects of changing various 
assumptions, including the basic assumption itself, on the results; we see how the results 
change if we neglect to subtract the background/foreground emission from the data, the 
effects of the emission of the dust associated with the H I gas, the effects for different values 
of the spectral emissivity index, the effects of a spatially varying 13 CO abundance, the effects 
of varying the column density per velocity interval or density, and how the signal-to-noise 
filtering has affected the results. 

3.1. The Effects of No Background/Foreground Subtraction 

As stated in Section 2.1 of Paper I, there are uncertainties in the subtraction of the 
large-scale emission (i.e., on the scale of the entire map shown in Figure 2 of Paper I). This 
large-scale emission was subtracted from the 140 fim, 240 fim, and H I maps. Even though 
this uncertainty was estimated to be 10%, it is still a good idea to see how this subtraction 
affects the model results. This was done by repeating the model in the LTE and LVG, 
one-component cases and in the LVG, two-component case (entire subsample) for the data 
without the background/foreground subtraction. The results were roughly similar to scaling 
up the data or, equivalently, scaling down the model curve. As such, the model results in 
the tested cases were roughly equivalent to those obtained for data that did indeed have 
the subtraction of the large-scale emission with a scale factor of about 0.9 applied to the 
model curve. Consequently, the systematic uncertainty in the observed data is much smaller 
than that of the total adopted calibration uncertainty of 40%. This therefore implies that any 
reasonable estimate of the uncertainty in determining the appropriate level of the large-scale 
emission to be subtracted will have an even smaller effect on the model results (equivalent 
to adjusting the scale factor by a few percent). 

3.2. The Effects of Dust Associated with HI 

Even though the HI has only small column densities on all lines of sight in the Orion 
fields (i.e. the average N(H I) is 5 x 10 20 cm -2 for those positions greater than 5-<r in 
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1^(140 yum), I„(240/xm), and I( 13 CO)), the dust associated with the HI may still have a 
non- negligible effect on the model results. To test this, the 140 /im/240 /im color tempera- 
ture of the dust was plotted against the atomic gas to molecular gas column density ratio, 
N(HI)/2N(H 2 ). Figure [JJ shows that the dust color temperature tends to about 21 K as the 
atomic to molecular gas ratio increases. The correlation between the color temperature and 
the atomic-to-molecular gas ratio has a confidence level of better than 99.99% according 
to the Spearman rank-order correlation test. (More specifically, the significance of the null 
hypothesis of zero correlation is less than 10~ 24 .) The curves represent the hypothetical case 
of having all the dust in the molecular gas at one fixed temperature and all the dust in the 
atomic gas at some other fixed temperature. The lower curve assumes that the dust asso- 
ciated with the molecular gas has a temperature of 16.5 K, while the upper curve assumes 
a dust temperature of 27 K for the molecular-gas-associated dust. Both curves assume that 
the atomic-gas-associated dust has a temperature of 22.5 K. Both curves together crudely 
describe the trends in the data. Consequently, both curves together imply that each line of 
sight either has molecular gas with cold dust, with Td = 16.5 K, or with warmer dust, with 
Tj = 27 K, along with atomic gas that has dust with a constant temperature of 22.5 K for 
every line of sight. The variations in dust color temperature would then be largely due to the 
variation in the atomic-to-molecular gas ratio (along with some scatter). This contradicts 
the picture represented by the models applied to the r 240 versus Td c plot. In that picture, the 
dust in the molecular gas does indeed vary in temperature from one line of sight to another, 
at least for the dominant component (i.e., component 1). 

To resolve this discrepancy, the simulated maps discussed previously were modified by 
adding a layer of HI and its associated dust with uniform properties throughout: a constant 
column density of 5xl0 20 H atoms cm -2 and a constant dust temperature of 22.5 K. Noise 
was added to the H I column density map that was the same as the value for the observed H I 
map. The 140 /im and 240 /im intensities for this H I layer were computed and added to the 
original noise-free maps. The noise for the new continuum maps, that include the H I layer's 
dust emission, was then recomputed from the prescription used previously (i.e., expressions 1 
and 2 of Paper II). The results are plotted in Figure [2] in the form of 140 /im/240 /im color 
temperature versus atomic-to-molecular gas ratio, analogous to Figure [JJ The simulated 
data does indeed reproduce the overall shape of the observed data, despite the lack of scat- 
ter in the former compared to the latter. This scatter in the simulated data can be increased 
realistically by including variations in the HI column density map and in its dust temper- 
ature. Nevertheless, it is clear that having only two possible constant dust temperatures 
in the dust associated with the molecular gas is not necessary for explaining the trends in 
Figure [TJ Figure [2] shows the same overall trends even though it uses model results that 
allow a changing dust temperature in the molecular gas. 
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The question remains as to the size of the change in the model results because of this 
HI layer and its dust. This question has already been answered indirectly in the previ- 
ous subsection, the subsection that dealt with the effects of no subtraction of the back- 
ground/foreground (i.e., large-scale) emission. Such a subtraction was necessary for the 
continuum maps (i.e., maps of dust emission) and for the HI map (i.e., map of atomic gas 
emission), but not for the 13 CO map. Therefore, the dust emission on the large scale is asso- 
ciated almost completely with the H I gas and not with the molecular gas. Accordingly, not 
subtracting the large-scale emission from the continuum maps is equivalent to piling on the 
atomic gas and its associated dust. In fact, it is equivalent to increasing the quantity of the 
H I- associated dust by factors of 4 to 5 over that in the background/foreground-subtracted 
maps. And yet, as described previously, the effect of this extra dust {not associated with 
molecular gas) was to change the results of the models (that only considered the dust as- 
sociated with molecular gas) in a way consistent with changing the scale factor by only 
about 10%. Therefore, the small amount of H I gas and its associated dust that remains in 
the background/foreground-subtracted maps will have an even smaller effect on the model 
results — equivalent to changing the scale factor by about 2 to 2.5%. It follows that the 
changes to the model results will be equivalent to only a tiny fraction of the full range of 
values for each parameter that is seen in such figures as Figure 21 of Paper I or in Table [2j 
(Or in the case where there are only lower limits, such as those for the densities listed in 
Table [21 those lower limits would be essentially unchanged.) 

3.2.1. The HI- Associated Dust and the Lower Temperature Limit of the Cold Dust and 

Gas 

One important point that remains is whether the dust associated with the H I is respon- 
sible for the model result that there is cold dust and gas at temperatures as low as 3 K. If 
so, then this result is incorrect. As discussed before, the H I-associated dust has little overall 
effect on the majority of the plotted points. However, this does not mean that the effect of 
this dust is negligible in every sub-grouping of points. Specifically, the vertical section of the 
model curve in the r 240 versus T dc plot (see Figures 20 and 24 of Paper I) that represents this 
very cold material is located between r 240 ~ 25 and about 50 MJy • sr -1 • (K • km • s -1 ) -1 for 
Tdc — 18 K. If the r 240 values were lower in this part of the plot, then no point would corre- 
spond to this "cold" section of the model curve. And the r 240 values would indeed be lower 
if the effect of this H I-associated dust were removed, thereby allowing a higher lower limit 
on the gas and dust temperature of component 1 — i.e., lower limits on T K1 and Tdi higher 
than 3K. Figure [3] shows plots of r 240 versus the HI fraction (i.e. N(H I)/[N(H I)+2N(H 2 )]) 
for the positions with Td c between 17 and 19.5 K and with signal-to-noise ratios greater than 
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or equal to 5 at 140 /im, 240 /xm, in the 13 CO J = 1 —>■ line, and greater than or equal to 
3 in the HI 21-cm line. These plots show a clear correlation (better than 99% confidence 
from the Spearman rank-order correlation test), suggesting that the higher r 240 values are 
indeed due to dust associated with the atomic hydrogen. The question is how strongly the 
H I-associated dust is contributing to the total 240 //m emission. This allows us to correct 
for the emission of this dust, thereby effectively giving us only the continuum emission and 
line emission from the molecular gas alone. 

Estimating the appropriate correction to r 240 for the HI associated dust is far from 
straightforward. One way is to make the crude assumption that the dust-emissivity per 
H-nucleon of the H I-associated dust, e HI , and for the H 2 -associated dust, e H2 , are each 
constant for all the positions in the subsample of points to be tested. This subsample is part 
of the sample used throughout this paper — signal-to-noise greater than or equal to 5 for 
L,(140/im), 1^(240 /xm), and I( 13 CO) — with the additional criteria mentioned above: 

a) Integrated intensity of the H I line greater than or equal to 3a. 

b) Tdc between 17 and 19.5 K. 

c) r 240 between 25 and 50 MJy • sr _1 • (K • km ■ s -1 ) -1 . 

The 240 /xm specific intensity, 1^(240 /xm), can be represented as 



Simple linear regression in three dimensions can be used to solve for the emissivities, e HI and 
e H2 , where we fit an equation of the form 



solving for the optimal values of the coefficients a and b. Since normal linear regression 
only uses the uncertainties in the z values, the fit should be repeated after interchanging 
the z values with the x values (and thereby using the ^-uncertainties) and then again after 
interchanging the original z values with the y values (and thereby using the ^/-uncertainties). 
This gives us three fits. For the first fit, the L,(240 /xm) are the z and the uncertainties in 
L,(240 /xm) are used in the fitting, N(HI) is x, and 2N(H 2 ) is y. The coefficients a and b 
then directly correspond to e HI and e H2 , respectively. The second fit has N(HI) as z, the 
uncertainties in N(HI) are used in the fitting, 1^(240 /xm) is x, and 2N(H 2 ) is still y. After 
solving for the optimum a and b values, expression (T5]) is rearranged to the form of ([1]) in 
order to solve for e HI and e H2 in terms of a and b. The third fit is similar to the second fit, but 
with 2N(H 2 ) as z, N(H I) as x, and I„(240 /xm) as y. Again, with the resultant a and b values, 



I„(240/xm) = e HI N(HI) + 2e H2 N(H 2 ) 



(1) 



z = ax + by 



(2) 
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the expression (j2J) is rearranged to the form of ([I]) and e HI and e H2 are found. (Another way 
to include the errors from all three quantities is to use the orthogonal regression method 
described at the beginning of Section 3 of Paper I.) In principle, the three different values 
determined for each of the emissivities can be used as a measure of their uncertainties. 

In practice, none of the fits were very good. The fits were applied using the molecular 
gas column densities (i.e., N(H 2 )) derived from the one- component, non-LTE models and 
those using those column densities derived from the two-component, two-subsample, non- 
LTE models. Typical reduced chi-square values were from 30 to 200, although one was 
as low as xl — 6. None of the fits using the two-component model column densities had 
reduced chi-square values less than 80. The poor quality of the fits reflects the invalidity of 
the assumption of spatially constant emissivities. If we nonetheless use the one fit with the 
least unacceptable xl value of about 6, then e H2 ~ 5e HI . Given that the maximum atomic 
gas fraction is about 0.3 (see Figure [3]), the H I-associated dust contributes less than about 
7% to the total 240 /im emission. If we tighten criterion c) above to r 240 between 40 and 
50MJy • sr _1 • (K • km • s -1 ) -1 to see how the results are changed, and if we ignore the fits 
with 10 or with a negative value of either e, then we find that e H2 ~ 9e HI . The least 
unacceptable fits suggest that the HI associated dust makes a contribution of about 3 to 
7% to the 240 /im emission. This result is similar, at least qualitatively, to the result found 
previously that the H I-associated dust emission has a negligible effect on whole sample of 
high signal-to-noise points. 

Given that no fit was acceptable, at least a simplistic correction should be applied to see 
what changes to the lower limit of the dust and gas temperature of component 1 are possible. 
The simplest kind of correction is to assume that the emissivities per H-nucleon are equal for 
the H I-associated dust and the H 2 -associated dust, a much more extreme correction than 
the fits to equation ([T]) would suggest. Accordingly, the observed I„(240/zm) values, and 
therefore the r 240 ratio values, must be scaled by 2N(H 2 )/[N(HI)+2N(H2)]. These correction 
factors were applied, using the one-component values for the molecular gas column densities. 
This is obviously not consistent with the two-component model results that predict the cold 
gas and dust. However, assuming correction factors from the one-component model column 
densities results in more extreme (i.e. further from unity) values of the correction factors 
that permit a greater appreciation of the tight constraint on the lower temperature limit. 
These correction factors were applied to the whole sample of high signal-to-noise positions - 
specifically, all the points with 1^(140 /im), I„(240 /im), and I( 13 CO) greater than 5a. Since 
the HI data were used, criterion a) reduced the sample to 609 points. The uncertainties 
in the correction factors were not propagated to the error bars in the r 240 values in the 
sample. The larger error bars in the r 240 ratios would have produced an ambiguity in the 
interpretation of the results: was the change in the lower temperature limit of component 1 
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really due to the correction for the H I-associated dust or was this simply due to the larger 
error bars? The results of fitting the two- component models to the corrected data should 
not be taken literally in any case and only point to the potential effects of correcting for the 
unwanted continuum emission. We are testing the hypothetical case of observing imaginary 
molecular clouds absolutely free of atomic hydrogen and its associated dust. 

The results of fitting a simple two-component model to the H I-corrected data are consis- 
tent with the previous two-component model results in Section 3.3 of Paper I. In particular, 
the results of these corrected models are within the range of values depicted in Figure 21 
of Paper I (ignoring the c a versus scale factor and the versus scale factor plots in fa- 
vor of the c Q ^- product versus scale factor plot), except that the density n cQ is as small 
as 10 cm" 3 . One noteworthy difference is that xl is higher for the two-component model 
fitted to the corrected data, xl = 6.41, than for the model fitted to the uncorrected data, 
xl = 5.69. If we now increase the lower temperature limit of component 1 from 2.8 K to 5 K 
and search for the optimum parameter values again, then xl increases to 7.54. Specifically, 
we set the lower limits of T K1 to 5 K and of AT to K, resulting in both T K1 and T^i > 5 K. 
Again keeping in mind that the points in the sample are not completely independent, the 
effective number of degrees of freedom is about 60. Adopting this number, the F-test tells 
us that raising the lower limit of component-1 temperature to 5 K can only be rejected at 
a confidence level of nearly 75%. Such a confidence level does not inspire much confidence, 
and we should probably accept the fit. Consequently, it is possible that gas and dust with 
temperatures between 3 and 5 K are not necessary for explaining the observations. If this 
is the case, then the correction factor of 1.6 to the one-component masses would have to be 
corrected downward. 

In general, it seems that the H I-associated dust cannot provide an alternative explana- 
tion for the cold dust and gas (i.e. temperature between 3 and 10 K). Nonetheless, there is a 
possibility that the lower limit of the component-1 temperature is about 5K instead of 3K. 
Then the estimated fraction of the total gas mass in the cold portion would be less than the 
original estimate of about 40-50%. 



3.3. Varying X( 13 CO) 



War in et al.l (119961 ) examined the photodissocation and rotational excitation of CO 



and its isotopologues in diffuse, translucent, and dense dark clouds. They found that the 
13 CO/ 12 CO abundance ratio (i.e. X( 13 CO)/X( 12 CO)) can vary by factors of 2 or 3 as a func- 
tion of depth into the cloud (see the left-most panel of their Figure 16). Given that T^ varies 
with N(H I + 2H 2 ) in the Orion clouds (e.g., see Fig 4 of Paper II), is it then possible that the 
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observed variation of r 240 with Td is because X( 13 CO) is varying with N(HI + 2H2)? Examin- 
ing the I( 13 CO) /I( 12 CO) as a function of N(H I + 2H 2 ) (an a pproximation of the ir Figure 16) 
shows no consistent trend and is at odds with the results of IWarin et al.l (119961 ). Therefore, 
it is very unlikely that variations of X( 13 CO)/X( 12 CO) could account for the overall trends 
seen in r 240 versus Td. 

Nevertheless, variations in X( 13 CO)/X( 12 CO) could still explain unusually small or large 
r 240 values for some points in the r 240 versus T d plots. In particular, the points near the top 
of the triangular cluster could be explained by a lower X( 13 CO). This is discussed in detail 
in Section H~5l 



3.4. Models with (3 ^ 2.0 



Johnstone fc Ballyl (119991 ) find that the dust emissivity index, (3, might have a range as 
extreme as from 1.5 to 2.5 in Orion A. Two-component models using adopted values of 1.5 and 
2.5 for (3 were applied to the data to see the effects on the results. The resultant parameter 
values were consistent with the ranges of values listed in Table [2], with one exception: the 
range of Tdo values is changed because (3 is different from before — all the inferred dust 
temperatures are changed. For (3 = 1.5, T d o = 21 K was found and, for (3 = 2.5, T d o = 16 K 
was the best fit. 

In short, with the exception noted above, the resultant parameter values are still within 
the range of values expected from the scale factor variations. Thus the most extreme varia- 
tions of (3 still largely give model results within the ranges listed in Table [2J 



3.5. Models with Varying N( 13 CO)/Av or Varying n(H 2 ) 

These were essentially one-component models with only three parameters on any one 
given line of sight: AT, N( 13 CO)/Av and n(H2). The approach here was to vary not only 
the Td and T K from one line of sight to another (keeping AT constant), but also to vary one 
of two parameters: the column density per velocity interval (N( 13 CO)/Av) or the density 
(n(H 2 )). To characterize all the lines of sight represented in the r 240 versus T d plot, there were 
9 parameters in total for a curve that was strongly decreasing through the triangular cluster 
of points at low T d , gently rising for the intermediate T d , and then strongly rising for high 
Td- Thus there were 3 intervals of T d that were delimited with 4 parameters: T 1; T 2 , T 3 , and 
T 4 . If we consider first only the models with varying N( 13 CO)/Av, then on the first interval, 
[Ti, T 2 ], the N( 13 CO)/Av was equal to jVj at 7\ and to N 2 at T 2 and the N( 13 CO)/Av value 
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at any point in this interval was determined by logarithmically interpolating between Ni 
and N 2 . The same approach was used on the last interval, [T 3 ,T 4 ], with N( 13 CO)/Av values 
equal to N 2 at T 3 and to N 3 at T 4 . On the interval [T 2 ,T 3 ], N( 13 CO)/Av was held constant 
at a value of N 2 . The fitted parameters were then AT, 7\, T 2 , T 3 , T 4 , Ni, N 2 , iV 3 , and the 
fixed density, n (with best-fit values of 0±1K, 16±1K, 17±1K, 24±1K, 27±1K, (5± J) x 
10 16 13 CO molecules ■ cm~ 2 ■ (km - s -1 ) -1 , (2± \) x 10 14 13 CO molecules ■ cm -2 • (km - s -1 ) -1 , 
(8 ± 5) x 10 16 13 COmofec«les • cm" 2 • (km • s" 1 ) -1 , and (3 ± \) x 10 3 H 2 molecules • cm" 3 , 
respectively). Achieving a reasonable- looking "best" fit presented a number of problems: 

1. The fitting is heavily biased towards the low-T d points (i.e., T d <21K). The points 
with T d >22K are almost excluded in the model fits. 

2. Even within the low-T d points, the fits tend to bypass the central area of the large 
triangular cluster of points. 

3. There is a very large hump in the model curve for T d around 22 K that overestimates 
the data points by more than an order of magnitude. 

4. The best fit for N 2 is more than order of magnitude lower than the observed large-scale 
N( 13 CO)/Av value (see Section 3.2 of Paper I). 

5. In order to fix the bias towards the low-T d points and remove the hump, it was necessary 
to apply weights to the data to reduce or remove this bias. This resulted in a chi-square 
value that was somewhat subjective, given that the choice of weights was somewhat 
subjective. 

The models with varying n(H2) had similar, though not identical, problems. For these 
models, the N( 13 CO)/Av parameters interchange with the n(H 2 ) parameters: AT, T 1; T 2 , 
T 3 , T 4 , ni, n 2 , n 3 , and the fixed column density per velocity interval, iV (with best-fit values 
of -2 ± 1 K, 14 ± 1 K, 18 ± 1 K, 22 ± 1 K, 30 ± 1 K, 10 ± * H 2 molecules ■ cm~ 3 , (5.6 ± \\ x 
10 4 H 2 molecules ■ cm~ 3 , 10 ± ® H 2 molecules ■ cm" " 3 , and (2 ± }) x 10 15 13 CO molecules ■ 
cm~ 2 • (km • s -1 ) -1 , respectively). As with the models mentioned in the previous paragraph, 
these models have biases that inhibit good fits through the T d >22K points or through the 
central part of the triangular cluster. However, the problems for n(H 2 )-varying models differ 
in three ways from those of the N( 13 CO)/Av- varying models: 

• The hump for the N( 13 CO)/Av- varying models centered at about T d = 22 K is much 
smaller for the n(H 2 )-varying models. This hump for the latter models only overesti- 
mates the data by about a factor of 2. 
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• The best-fit N( 13 CO)/Av value is still low, but is now within a factor of 2 of the lower 
limit imposed by the large-scale observed N( 13 CO)/Av value. 

• The best way to reduce the bias towards the low-T d points was to fix the Ti, T 2 , T 3 , and 
T 4 values. The choice of these values and their uncertainties was somewhat subjective. 

The best-fitting model curve gave a reduced chi-square of 9.2. 

While the best of these models is much better than the one-component models that 
used the basic assumption (i.e., xl — 16.5 and 16.9), it is still much worse than the two- 
component models (i.e., xl = 5-7 and 5.3). In addition, the models with varying n(H2) (or 
varying N( 13 CO)/Av) require subjective judgement in determining the values to fix for some 
of the parameters (or in adjusting the weights of different subsamples of data points). And de- 
spite having fewer parameters (9 as opposed to 15) than the two- component, two-subsample 
models — the most elaborate of the models that obey the basic assumption (within each sub- 
sample at least) — the models with varying n(H 2 ) or varying N( 13 CO)/Av are actually more 
complicated than these two-component models. They are more complicated because these 
models essentially divide up the sample of data points into three subsamples instead of two, 
and two of these subsamples (in intervals [T 1 ,T 2 ] and [T 3 ,T 4 ]) allow n(H 2 ) (or N( 13 CO)/Av) 
to vary in addition to Tj (along with T K so as to keep AT constant). In contrast, the 
two-component, two-subsample models have only Td varying (along with T K so as to keep 
AT constant) within each subsample — the basic assumption is obeyed in each subsample. 

In short, relaxing the basic assumption by using models with simple systematic varia- 
tions of either N( 13 CO)/Av or of n(H 2 ) results in poorer, and more subjective, fits than the 
best models following this assumption. 

3.6. Signal-to-Noise Considerations 

Given that the sample of points modeled in the r 240 versus T dc plots represents only 
25.8% of the Orion Fields, how well does this sample represent the physical conditions in 
the Orion clouds as a whole? A related question is what fraction of the positions that have 
gas/dust emission from the Orion clouds is represented by the sample? The contour maps 
in Figure 2 of Paper I suggest that roughly three-quarters of the Orion Fields have gas and 
dust. Some low- level emission could come from elsewhere along the line of sight, so maybe 
only half of the Orion Fields are occupied by the Orion clouds. However, even in this extreme 
possibility, the modeled positions still only represent about 50% of the Orion clouds' area 
- still not a majority of the positions. Clearly, just how strongly the inferred physical 
conditions depend on the signal-to-noise threshold used in sample selection must be tested. 
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The most straightforward test, namely lowering the signal-to- noise threshold until most 
of the Orion Fields are represented, is not practical. There would be 2 to 4 times the 
number of points to fit and the majority of these would have huge error bars. Consequently, 
this sample of points would be difficult to model reliably. The alternative test is to raise 
the signal-to-noise threshold and see how the inferred physical conditions change with this 
threshold. This has the disadvantage of going from a minority of the points to an even 
smaller minority. Nonetheless, this alternative has the very strong advantage that only very 
high signal-to-noise points are modeled, thereby yielding more reliable fit parameters. 

To do this alternative test, the threshold was increased appreciably: from 5a to 20a in 
Ij/(140 /mi), Ijy(240 /mi), I( 13 CO). This sample has only 6.6% of the Orion Field positions. An 
LVG, two-component model was fitted to this sample. The resultant parameter values were 
consistent with those specified in Tabled with one notable exception: T d0 was 17 K instead 
of the usual 18 K. This change in the component-0 dust temperature is not surprising. The 
signal-to-noise ratio is proportional to the surface brightness, and this depends mostly on 
the gas/dust column density. So signal-to- noise ratio is roughly equivalent to depth into 
the clouds. Increasing this ratio's threshold is almost like filtering out the cloud edges and 
looking more deeply into the clouds. The T d0 is the roughly constant temperature on the 
scale of the Orion clouds of some dust component. If this dust component is on the surfaces 
of the clouds, then it is the temperature of the dust heated primarily by the general ISRF. 
If this dust component is just below the clouds' surfaces, then this is the temperature of the 
dust heated primarily by an ISRF shielded by the surface layers of gas and dust, resulting in a 
lower temperature. Then by extrapolation, modeling of all the positions in the Orion Fields 
(if there were sufficient signal-to-noise) could yield T d0 — 19 K. 

Accordingly, modeling the full spread of points in the r 240 versus T dc plot probably 
requires a spread in the physical parameter values that have been held constant. For example, 
the horizontal spread of the triangular cluster in an r 240 versus T dc plot (e.g., see Figure 20 
of Paper I) probably means a component-0 dust temperature varying between about 16 and 
19 K. (As such, an updated version of Table 2 of Paper II becomes Table [2] of the current 
paper, which includes this estimate of the range of T d0 values.) Analogously, other parameter 
values, such as the densities and column densities per velocity interval, must also vary to 
"fill" the space occupied by the sample of points, as suggested in Section 3.5 of Paper I. 
Nonetheless, the current modeling effort is a sufficient first approximation. 

The most important result of varying the signal-to-noise ratio is that the sample with 
only 25.8% of the Orion Field positions may indeed represent the bulk of the Orion clouds. 
There was no appreciable change (except a slight change in T d o) in going from the factor 
of ~4 from 6.6% of the fields to 25.8%. Therefore, extrapolating the extra factor of ~4 to 
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100% may also yield no appreciable change. This is far from certain, of course, because 
the extrapolation from 25.8% to 100% could cross some depth threshold that appreciably 
alters the parameter values. On the other hand, the effective spatial resolution of these 
observations is 8 pc at the distance of the Orion clouds. Averaging over such large size scales 
may reduce the effects of varying the sample. 

In short, varying the signal-to-ratio used in selecting the sample may have little effect 
on the derived gas/dust physical conditions. 



4. Scientific Implications and Discussion 

The astronomical implications of the results are numerous and include modifications 
to models of dust /gas thermal coupling, to estimates of mass and kinetic temperature of 
molecular gas on galactic scales, and to our understanding of the X-factor to name a few. 
These issues and others will be discussed in the following sections (or in subsequent papers 
as for the X-factor), after discussing the appropriateness of using the LVG models. 



4.1. Use of the LVG Code 

Some literature suggests that the LVG models can give inconsiste nt results and tha t pho 



todissociation region (PDR) models remove the inconsistencies (e.g.. lMao et al.l 120001 ) . The 
PDR models represent a nearly complete explanation for the emission strength of molecular 
lines, whereas the LVG models merely relate the molecular line strength to simple physical 
parameters in a simplified case. Consequently, the PDR models should yield more reliable 
estimates of physical parameters like temperature and density than those of the LVG mod- 
els. However, some papers that use CO line ratios to claim that the LVG model results are 
unsatisfactory compared with those of PDR models often suffer from flawed and inconsistent 
arguments. 



Mao et al.l (120001 ). for example, model the physical conditions in the central ~500pc 
of the galaxy M82. They claim that their CO data yield physically unreasonable results 
when using the LVG models — results that supposedly become more reasonable when using 
PDR models. The problems with this claim is summarized below. They state, among other 
things, that the LVG models imply gas densities that are too low, cloud sizes that are too 
large, and area filling factors that are inconsistent with volume filling factors. Their first 
claim of low gas densities is connected with their restriction of the X( 13 CO)/(dv/dr) value. 
They use the observed large-scale velocity width of the lines compared with the size of 
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the observed region to fix the numerical value of X( 13 CO)/(dv/dr) within a narrow range. 
However, the large-scale (i.e., on the scale of the beam) X( 13 CO)/(dv/dr) has little to do 
with this parameter value within the individual clumps responsible for much of the observed 
emission. If their approach were also applied to the volume density, then it would have been 
equivalent to dividing the mass of gas within the beam by the volume of this gas to estimate 
and fix the gas density, which is well known to be only a rough lower limit to the density 
within the clumps. They did not make this mistake, but did make the equivalent mistake for 
the X( 13 CO)/(dv/dr) value. Had they allowed X( 13 CO)/(dv/dr) to vary over a wider range, 
they would have satisfied the density constraints suggested by the other observations they 
mentioned. Their second claim of overly large cloud sizes implied by the LVG model results 
was based on two different methods. Both methods, however, overestimated the cloud sizes 
by the ratio of the observed line velocity width to the cloud velocity width, which is about 
an order of magnitude. One method was dividing the column density of a single cloud by 
its volume density, where both the column density and volume density were LVG results. 
The second method estimated volume and area filling factors from LVG results, their ratio 
giving the size of a single cloud. Both methods were applied incorrectly. The first method 
converted from N( 13 CO)/Av to cloud N(H 2 ) by using the entire observed line width, which 
also includes the rotation of the observed galaxy, instead of using an estimate of a cloud line 
width (i.e. 5-10 km ■ s^ 1 ). The second method suffered from the same overestimate. The 
cloud size depends on the ratio of the volume to area filling factors, but the area filling factor 
must be for all the gas at all the velocities within the line profile (assuming that there is 
never more than one cloud on any line of sight for all the velocities within the line profile). 
However, they clearly used the ratio of the observed line strengths to those of the model, 
which gives the area filling factor within a narrow velocity interval and not the area filling 
factor over the entire line profile. This second method must include the ratio of the cloud 
velocity width to the line velocity width, again reducing the estimated cloud size by an order 
of magnitude. The third claim of inconsistent area and volume filling factors is weak at best, 
given that the relationship between the two is not a fixed, straight forward expressio n that 
applies in every case. If we nonetheless accept the expression used by lMao et al.1 (120001 ) . then 
there was indeed a minor dis crepancy between the two types of filling factor for the LVG 
model results. However, what iMao et al.1 (120001 ) ignored entirely was that the corresponding 
discrepancy for the PDR models is much larger than that for the LVG model results. The 
filling factor argument was applied in a clearly biased manner. In short, all three claims are 
based on arguments that are faulty or biased or both. 



In addition to the problems above, t here is the s t rong evidence provided by lWeiss et al. 



( 2001 ): they used data similar to those of Mao et al. ( 200ol) and recovered th e sa me physical 



conditions, while using only LVG models. Obviously, if I Weiss et al.l (120011 ) and IMao et al. 
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( 120001 ) agree on the physical conditions, then they dis agree on the neces sity of the PDR 
models for recovering those conditions. And given that IWeiss et al.l ( 20011) recovered those 
conditions with only the LVG models, then the claim of iMao et al.l (120001 ) that PDR models 
are necessary is clearly incorrect. Therefore, the LVG models are clearly as reliable as the 
PDR models when using CO lines (at least on scales of a few or more parsecs). 

Consequently, the LVG models applied to the large-scale physical conditions in the Orion 
molecular clouds, as described in this paper, are also as reliable as the PDR models. 



4.2. Comparison of Derived T K , n(H 2 ), and N( 13 CO)/Av with Previous Work 



There are few papers that discuss the molecular gas physical conditions of the entire 
Orion A and B clouds, as inferred from two or more transitions of CO. ISakamoto et al.l (119941 ) 
inferred these physical conditions using the J = 2 — > 1 and J = 1 — > lines of 12 CO. Of 
all the rotational lines of 12 CO and its isotopologues, these are the least sensitive to the gas 
physical conditions. Nevertheless, they at least provide a rough comparison with the physical 
conditions obtained in the current paper that used the dust-continuum to gas-line ratio. As 
discussed in Paper II, the 13 CO J = 1 —>■ emission is dominated by that of component 1 
in all but the few points with Tj = T K ^4K in that component. Since this is molecular gas 
emission, this emission, rather than the FIR continuum, better identi fies the component wit h 
which to compare single-component model results, such as those of ISakamoto et al . (1994). 



Conse quently, the component-1 parameters in Table[2]are compared with the lSakamoto et al. 
rtl994T ) results. 



The physical conditions inferred by ISakamoto et al.l (119941 ) are consistent with those 
inferred in the current paper. They found the following physical conditions in the Orion 
molecular clouds: 



N( 12 CO)/Av between about 1 x 10 16 and about 3 x 10 18 12 CO molecules ■ cm 2 , which 
corresponds to N( 13 CO)/Av between about 2 x 10 14 and 5 x 10 16 13 CO molecules ■ cm -2 , 

n(H2)<^3 x 10 3 cm~ 3 over much of the clouds' areas, except for n(H2) - 2 x 10 2 cm -3 
in the clouds' peripheries, 

T K between 10 K and 40 K, the latter temperature found near the HII regions. 



Sakamoto et al.l (119941 ) chose to do their LVG analysis with diagrams of constant 
X( 12 CO)/(dv/dr) instead of constant N( 12 CO)/Av. They considered only values of 1 x 10~ 4 



-20- 



and 1 x 10 5 (km • s 1 • pc *) 1 . Given that the 13 CO J = 1 — > cl umps have veloc- 
ity w idths and sizes consistent with dv/dr of a fewkm ■ s _1 • pc -1 (see iNagahama et al. 



1 • pc- 1 ; 



19981 ). the maximum X( 12 CO)/(dv/dr) value should be a few x 10~ 5 (km ■ s 
given the 12 CO abundance mentioned in Appendix A of Paper I. Also, given that some 
structures ha ve dv/dr of around 10 km • s _1 • pc -1 — even though seen with other tracers 
like CH3OH (ICernicharo et al.lll999l ) these structures are nonetheless real — the minimum 
X( 12 CO)/(dv/dr) should be around 10 -6 (km • s _1 • pc -1 ) -1 . Nevertheless, because such 
structures have higher densities than those normally inferred from CO obse rvations on par- 



sec sca les, the range of N( CO) / Av values are probably still roughly those of lSakamoto et al. 



(1 1994 ). This range is larger than the range of 



N c i( la CO) 



values listed i n Table El but includes 



the range listed in that table. The upper limits of N( 13 CO)/Av for ISakamoto et al.l (119941 ) 
and for Tableware similar to within a factor of 3, but the lowe r limits differ by more than 
an order of magnitude. This disagreement is probably because ISakamoto et al.l (119941 ) did 
not use the large-scale N( 12 CO)/Av of the entire cloud as a rough lower limit on that of 
the clumps. Their densities are consistent with the lower limit in Table [2] for component 1, 
except for their density for the peripheries. This density is consistent with the lower limit we 
found for component 0. Accordingly, these peripheral regions may represent a low-density 
envelope surrounding the entire clouds and are seen in projection along the clouds' edges. 
The clumps in this envelope may be, in fact, the component of the current paper. 



The comparison between the kinetic temperatures of lSakamoto et al.l (119941 ) and those 
of the current paper is complicated by the high optical depths of the 12 CO lines observed by 
the former. The high optical depths imply that the inferred densities and column densities 
per velocity interval are biased towards of the surfaces of clumps. This is especially true for 
the inferred kinetic temperatures, because, in the optically thick case, these temperatures are 
more directly related to the observed line radiation temperatures than are the other physical 
parameters. The high optical depths then mean that the warmer component will dominate 
the 12 CO emission in these lines. We then simplistically assume that the T^c < 20 K subsam- 
ple is essentially component and that the Td c > 20 K subsample is component 1. (Note that 
this is not consistent with the choice of component used to compare the N( 13 CO)/Av and 
n(H 2 ) results. If these other physical parameters had been compared with the component-0 
results instead of those of component 1, then the agreement would still have been reasonable, 
given the loose restriction on the parameters.) With this approximation, the gas kinetic tem- 
peratures as sampled b y the 12 CO are approxim ately the same as the 140 /xm/240 /im dust 



color temperature, Td c - ISakamoto et al.l (119941 ) found T K = 40 K near the HII regions and 
10-20 K away from the H II regions. The current paper finds T^ = T K ~ 25 K near the H II 
regions an d 15-20 K away f r om th ose regions, a range of values less extreme than, and inside 
of, that of lSakamoto et al.l (I1994I ). The range is less extreme in the latter case because the 
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spatial resolution is worse by factors of 6 to 7 than those of ISakamoto et al.l (119941 ) 



In short, the physica l conditions derived here basically agree with the previous results 
of Sakamoto et ali ( 1994 ). 



4.3. Column Density Determinations and the Two-Component Models 

As found in Paper I (and Paper la), the column densities derived from the 13 CO J = 
1 — > line agreed significantly (at a confidence level better than 99.9%) better with those 
derived from the far-IR continuum for the two-component models than for the one-component 
models. An important difference between the two components is that component has 
unvarying gas and dust temperatures from sightline to sightline, whereas component 1 has 
these temperatures varying spatially (while maintaining a constant dust-gas temperature 
difference). Thus each sightline has two gas temperatures (and two dust temperatures), 
except for those few sightlines where the component-0 and component- 1 temperatures were 
the same. Hence, two temperatures along each sightline was necessary for recovering reliable 
column densities. 



This result is similar to that found by lSchnee et al.l (120061 ). They compared the continuum- 
derived column densities with those derived from extinction for the Perseus and Ophiuchus 
molecular clouds. The scatter in those plots more or less matched that in a simulated cloud 
that assumed isothermal dust on each sightline. They concluded, therefore, that deriving 
reliable column densities requires assuming variations of temperature along each sightline. 
This supports the result of the current work that at least two temperatures are needed on 
each sightline for estimating reliable column densities. 



4.4. AT = 

The result that the gas and dust temperatures are the same is unexpected, both theoret- 
ically and observationally. For example, t he theoretical model of PDR s applied to the Orion 



Nebula and its associated molecular gas (ITielens fc Hollenbachlll985l Jal) predicts AT ~ —70 



to +20 K for cloud depths for which the dominant form of carbon is CO (i.e. A v >3mag). 
This model, however, does not apply to multi-parsec scales in the Orion clouds: it uses a 
far-UV radiation field strength of G Q = 10 5 and a density of 2.3 x 10 5 H nuclei ■ cm' 3 ; both 
are too high for the molecular clouds on larger scales. The far-UV radiation field on such 
scales is roughly G Q ~ few (see Figure 17 and Section 14.11 of W96) and the density could 
be as low as few x 10 3 H2 cm -3 (see Table 1). Consequently, the "standard" PDR model 
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adopted by IMochizuki fc Nakagawal (120001 ) — which assumes n(H2) ~ 10 3 cm -3 , at cloud 
depths for which H 2 is the dominant form of hydrogen (i.e. N(H 2 )^ 10 21 cm -2 ), and G Q = 10 
- is more appropriate for comparison with the work done here and yields AT ~ —13 to 
+5 K. Considering the uncertainty in the AT of the current work, these theoretical val- 
ues are factors of about 5 to 10 too large. Given tha t there are many heating and cooling 
mechanisms in PDRs (see lTielens fc Hollenbachlll985l . for a comprehensive description) with 
different dependences on the density and on the radiation field (as well as on other quan- 
tities), adjusting the theoretical expressions for these many mechanisms to give AT near 
zero would only yield the desired result for an improbably narrow range of physical condi- 
tions (i.e., narrow range of G Q and n(H 2 ) values). On the other hand, simply increasing the 
gas-dust thermal coupling by the required fact or will easily achieve the de sired result. This 
corresponds to the A gd function developed by iBurke fc Hollenbach (|l983h and represents a 
cooling mechanism of the gas due to its collisions with the dust (and, of course if A gd < 0, 
it represents a heating mechanism of the gas due to those collisions). Therefore, increasing 
Agd by factors of 5 to 10 can explain the current observations. 



A commonly used form o f A g n is that of lGoldsmith (2001) and is g iven a s equation (lAlOO 
in Appendix [A] iGoldsmithl (120011 ) used the IBurke fc Hollenbach! (119831 ) expression and 
adopted certain parameter values, including a grain size that was a little too large. Also the 
Burke Hollenbach! (119831 ) expression only assumes a single grain size , and not the rang e 
of grain sizes that exists in the ISM (see for example iMathis et al.lll977l ; iDesert et al.l ll990'l. 
When a reasonably realistic range of grain sizes is considered, the A g d of Goldsmith (j2001 ) 
is increased by factors of 3 to 4 (see Appendix |AJ. This is not quite the factor of 5 to 
10 desired, but additional increases are easily possible when one considers grains with non- 
spherical shapes or with projections on their surfaces. As mentioned in Appendix [A] A g( j is 
proportional to the ratio of the grain geometric cross-section to the grain volume (assuming 
uniform grain density). If the grains are elongated, then it would be easy to increase this 
ratio. (In reality it is the cross-section averaged over all viewing angles that is important 
here. However, even modest elongations of a factor of a few would still result in a larger 
average cross-section than a sphere with the same volume.) Alternatively, projections on the 
grain surface could also increase this ratio, but, since the relevant area is the cross-section 
rather than the total surface area, these projections would have to be large compared to 
the grain size. In any event, achieving an additional factor of 2 is possible. This would 
mean that simple geometric considerations could increase the commonly used form of A g d 
by factors of 6 to 8. Therefore, A g d can indeed be larger than had been previously assumed 
and could possibly explain the AT ~ result for the Orion clouds. 



Observationally , AT ~ is unexpected as well (e.g., IWu fc Evanslll989l ; iMangum et al. 



19991 ; iLis et al.ll200ll ). As discussed in the introduction, the different temperature and density 
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sensitivities of the dust continuum emission and gas line emission can result in incorrect 
inferences of the relative dust and gas temperatures at each point along the line of sight; 
the continuum and line emission preferentially trace different regions of the ISM within the 
same line of sight. In addition, in many cases the uncertainty in the dust temperature or 
in the gas tempera ture, or both, is large enough that Td = T K cannot be ruled out. In 



Wu fc Evand (119891 ). for example, the uncertainties in Td were usually ±1 K or ±2K. The 



uncertainties in T K were not explicitly listed for some sources, but would be at least about 
10% due to the stated calibration uncertainty. Even though they did observe two lines of 
CO — J = 1 — > and J = 2 —>■ 1 of 12 CO — they did not use their ratio to estimate T K ; as 
shown in Figure 1 in the introduction of Paper I, this gives only a very uninteresting lower 
limit on T K . Instead, they used the peak radiation temperature of each line, Rayleigh- Jeans 
corrected and corrected for the cosmic background, to estimate T K . They assumed the lines 
to be optically thick and thermalized. However, they also implicitly assumed that the gas 
fills the beam at the peaks of the lines. As we have found in the current paper, this is 
not necessari l y true , and may not be true even with the superior angular resolution of the 



Wu fc Evand (119891 ) observations (a factor of 30 to 60 smaller beam). In addition, they used 



the 60 fim and 100 /zm observations of IRAS. As stated before, the 60 /zm emission suffers the 



contamination of emission fro m stochastica l ly he ated dust grains (e.g., IDesert et al.lll990l . 
W96). However, because the IWu &: Evand (119891 ) observations are on angular scales of V 



to 2', the radiation fields on such small size scales could be large enough that the 60 fim 
emission largely comes from grains in thermal equilibrium (i.e., the grains that would be 
stochastically heated in a normal interstellar radiation field reach thermal equilibrium in a 
strong radiation field). A crude extrapolation of the trend in the data in Figure 10a of W96 
suggests that the column densities derived from 60 and 100 /im data will agree with those 
derive d from longer wav elengths for Td — 60 K. Since Td = 18 K for the general ISRF (see 
W96, IDesert et al.lll990l ). and given that the ISRF is proportional to T^ +4 , then radiation 
fields of at least G Q ~ few x 10 2 ensure that the 60 //m e mission originates la rgely from 
grains in thermal equilibrium. One of the sources observed by lWu fc Evand (Il989l ). B35, has 
a radiation field of G Q ~ 30 (jWolfire et al.lll989l . and references therein), roughly an order 
of magnitude too low to exclude the likelihood of stochastically heated grains contributing 
to the 60 /im emission. Therefore, for the source B35, and probably a few others in their 
list, the derived dust temperature overestimates the dust temperature of the large thermal 
equilibrium grains. In short, their estimates of T K are lower limits and their Td values are 
likely to be overestimates (especially in the case of B35). Accordingly, their conclusion that 
the observations are consistent with Td > Tk for the majority of their sources still does not 
exclude the possibility of Td = T K for these same sources. 



The same can be said for the observations of iMangum et al.l (119991 ). They observed lines 
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of formaldehyde, H2CO, towards dense gas condensations in NGC2024 in the Orion B cloud 
with angular resolutions of 12", 19", and 30". They found T K values from around 50 K to 
ar ound 250 K. They th en compared the derived gas T K values with the dust temperatures 
of iMezger et al.l (119921 ). derived from the ratio of the 870 /im- to the 1300 /im-continuum 
emission at resolutions comparable to those of the H2CO line observations (i.e., 24" and 
8", respectively). These continuum observations imply Td = 19 K. While, at face value, a 
difference between a T of 250 K and a Th of 19 K may seem substantial, the uncertainties 
in T K (see Table 4 of iMangum et al.l Il999f ) suggest that T K — T d is significant at levels of 
only 2 to 4cr. These would be satisfactory levels of significance, except that the correct 
uncertainty in T K — T<j m ust also inclu d e the uncertainty in Td as well — an uncertainty 
that was ignored entirely. IMezger et al.l (119921 ) state a 20% uncertainty in their continuum 
fluxes, implying an uncertainty of about 30% in the 870 /im to 1300 /im continuum. This 
uncertainty implies that Td = 19 K is consistent with Td = 7K to 00. (Note that, even if we 
optimistically assume the continuum intensity ratio uncertainty to be only 20%, the upper 
limit on Td would still be 00.) In other words, the continuum observations do not place any 
upper limit on T<j. In fact, the difference between 19 K and 250 K corresponds to only a 16% 
change in the 870 /im/1300 /im intensity ratio — a change of ~0.5-o\ Therefore, again, the 
observations do not exclude the possibility of Td = T K . 

In c o ntrast , observations of giant molecular cloud cores in the Galactic Center by 



Lis et al.l (120011 ) seem to genuinely rule out equal dust and gas temperatures. They use 
continuum observations at a number of wavelengths from 45 to 850 /im and observations of 
the molecular lines of H 2 CO, CS and other molecules to determine reliable dust and gas 
temperatures. They find two components of dust: a warm component with Td — 35 K, 
which dominates for A^ 100 /im, and a cooler component with Td — 18 K, which dominates 
for A^ 100 /im. (Note that they also estimated the radiation field strength to be G Q ~ 500 
to 1000. Therefore, the shorter wavelengths can also give reliable dust temperatures in this 
case.) The molecular gas temperatures are T K ~ 60 to 90 K, implying gas-dust tempera- 
ture difference as high as about 70 K. However, given the appreciable foreground emission 
towards the Galactic Center and the lack of velocity information from the continuum obser- 
vations, contamination of these continuum observations by such foreground emission cannot 
be ruled out. This foreground e mission would be from the dust in the Galactic disk and 
has a temperature of about 18 K (ISodroski et al.lll994j ). Consequently, only the 35 K dust 
might be directly associated with the observed cores in the Galactic Center. The molecules 
observed have transitions with high critical densities {n crU > 10 4 H 2 ■ cm -3 ) and it could be 
argued that the observed transitions are only sampling the densest portion of the molecular 
gas. However, gi ven that the bulk of the gas in the Galactic Center is high-density gas 
( IBally et al.lll987l ). the observed transitions are probably sampling most of the molecular 
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gas in the observed cloud core. Therefore, the gas-dust temperature difference cannot be 
any smaller than about 25 K, but is probably smaller than 70 K. 

Considering the observations, we can draw an important conclusion about AT: it seems 
that AT ~ is not excluded for a weak ISRF of G Q ^, 10 2 and that it is excluded for a 
strong ISRF. Another possibility is that, given that the Galactic Center represents a unique 
Galactic environment, finding AT ^ may have more to do with other physical conditions 
than simply the strength of the radiation field and the gas density. Nevertheless, for now, 
a good working assumption is that AT is indeed near for G Q < 10 2 and that AT is quite 
different from that, i.e. | AT|> 25 K, for G ^, 10 3 . If this assumption is correct, then AT ~ 
on multi-parsec scales for most molecular gas in the Galaxy, and AT ^ on these scales 
in the Galactic Center or in regions with large-scale star formation. This has a number of 
consequences, including the following: 



Galactic-scale molecular gas temperatures are nearly double the temperatures previ- 
ously believed. Applying corrections for the cosmic background and for the Rayleigh- 
Jeans approximation, the peak radiation tem peratures found for t he 12 CO J = 1 — > 
line in large-scale surveys of the Galaxy (e.g., ISanders et al.lll985l ) suggest that T K ~ 
10 K. If AT is indeed close to zero, then the tr ue T K is close to that of Td on large 
scales, which is T d ~ 20 K (ISodroski et al.lll994l ). As discussed in Section 3 of Paper I 
and illustrated in Figures 11, 18, 23, and 26 of Paper I, the molecular gas sampled by 
the 12 CO J = 1 —>■ line does not fill the beam within each velocity interval within the 
line profile, especially if the linear beam size at the source is parsecs. 



The gas not completely filling the beam in each velocity interval may better explain the 
X-factor. The usual explanatio n given for the N(H2)/I(CO) factor is some variation 
of that of iDickman et al.l (119861 ). that molecular clouds are virialized and that the line 
width indicates cloud mass and, therefore, cloud column density. In fact, in some cases 
the velocity- widths o f a cloud are only weakly correlated with column densities (e.g., 



sec 



Heyer et al.lll996l ). Consequently, a better explanation sprin gs from hav ing filling 



factors less than unity. This will be discussed in a future paper (jWallll2007bl ). 



3. AT = constrains proposed explanations of the dust grain alignment that has been ob 



served in the ISM flHiltnerlll949 



Hallll949h . For example, the Davis-Gree nstein mecha- 



nism is the relaxation of paramagnetic grains spinning in a mag netic field (jDavis fc Greenstein 



19511 ) . This relaxation mechanism requires that Tj ^ T K ( iJones Sz Spitzerl 119671 ) . 



However, there are number of other possible mechanisms that could ex plain dust grain 



alignment that require no such difference in temperatures (e.g., see lLazarian et al. 



19971 : lAbbas et al.ll2004l . and references therein). 
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Therefore, having equal gas and dust temperatures has a number of interesting conse- 
quences that are not necessarily contradicted by theory or observations. 



4.5. Cold Gas/Dust 

The two-component masses compared with those for the one-component masses in Ta- 
ble 6 of Paper I imply about 60% more mass of gas and dust in the Galaxy than previous 
estimates suggest. These es timates are on the order o f 5 x 1O 9 M of gas (i.e. molecular 



and atomic) in the Galaxy (jDame!ll993l : ISanderslll993l ). If the model results for Orion are 



taken at face value and if these results apply to other clouds throughout the Galaxy, then 
this total gas mass increases to about 8 x 10 9 M . This increase is due to some positions in 
component 1 having temperatures below 10 K and as low as about 3 K, nearly that of the 
cosmic background. Such cold dust and gas emits only weakly per unit mass, allowing much 
gas and dust to be "hidden" for the observed brightness. Indeed, Table 7 of Paper I lists the 
cold gas mass and it is about 40% of the total mass listed in Table 6 of Paper I (adopting 
Case 4 as the more realistic of the two listed in Table 7 of Paper I). Accordingly, the total 
gas mass is the warm gas mass increased by about 60% to allow for the cold gas mass. Such 
an increase, especially if it applies to the entire ISM of the Galaxy, is substantial and its 
validity must be examined carefully. 

As mentioned in Section 3.6 of Paper I, the existence of this cold dust and gas depends 
on the basic assumption used in the modeling. Even if this assumption has provided a good 
physical description for most of the points in the r 240 versus Td c plots (e.g., see Figures 20 
and 24 of Paper I), it does not necessarily apply to all of the triangular cluster of points 
from Td c — 15 to 21 K and r 240 ~ 10 to 70 MJy • sr -1 • (K - km - s -1 ) -1 , especially to those with 
r 240 > 30 MJy ■ sr" 1 • (K • km ■ a -1 ) -1 . Indeed, it is the r 240 > 30 MJy ■ sr" 1 ■ (K • km ■ s -1 ) -1 , 
Tdc = 18 K points that require Td as low as 3 K when the basic assumption applies. If these 
high-r 240 , T dc = 18 K points are explained by other means, then such cold gas and dust may 
not be present. In other words, let us abandon, for the moment, all models that use any form 
of the basic assumption, at least for these points. There then exist a number of possibilities 
for the high-r 240 , Td c = 18 K points: 

• The N( 13 CO)/Av or n(H2) is different from those of the rest of the points. Models with 
N( 13 CO)/Av or n(H 2 ) that vary smoothly with T dc were discussed in Subsection 13.41 
Even though such models have more difficulties than models using the basic assump- 
tion, it does not exclude this possibility. The N( 13 CO)/Av value would be higher or 
the n(H2) value would be lower. As we saw in Subsection I3.4[ n(H2) would be as low 
as about 10 cm -3 . Figure 29 of Paper I shows that these points occur mostly on some 
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edges of the Orion clouds. A lower density for these points is consistent with them be- 
ing on the cloud edges. However, such low densities imply a peak T R ( 13 CO J = 1 — > 0) 
either barely as strong as, or as much as an order of magnitude weaker than, the ob- 
served 13 CO J = 1 — > line strength. The other possibility of higher N( 13 CO)/Av 
would indeed give strong 13 CO J = 1 — > emission, but implies N( 13 CO)/Av on the 
cloud edges a factor of a few higher than that for the cloud central regions. This is 
possible because N( 13 CO)/Av is not equivalent to N(H 2 ), but still seems somewhat 
implausible. 

The points with high r 240 values at T dc ~ 18 K have appreciable emission of dust 
associated with atomic hydrogen. Given that this gas is largely found at some cloud 
edges, this explanation seems reasonable; the atomic hydrogen and its associated dust 
would be found on molecular cloud edges, thereby providing the shielding necessary for 
the existence of the molecular gas. However, we examined this in Subsection 13.21 and 
found, based on the current data, that such emission was not likely to be important. 
The H I-associated dust contributes negligibly to the observed r 240 values. Nonetheless, 
as discussed in that section, the gas and dust of component 1 with temperatures 
between 3 and 5K might still be explained by H I-associated dust. 

The AT is different for the high-r 240 , T d = 18 K points. Figure 7 of Paper I illustrates 
that AT varying smoothly from to less than — 16 K (e.g., ~ — 20 K) could account 
for the vertical extent of the points at Td = 18 K in the r 240 versus Td plot. However, 
Figure 17 of W96 suggests that G„ is only a few and this would give AT less extreme 



than those of iMochizuki fc Nakagawal (120001 ). rendering AT ~ —16 to — 20 K unlikely. 



Figure 7 of Paper I also suggests that AT = +14 K for these points could also account 
for their vertical extent. But again, this is too extreme for the given G Q . Also, this 
AT combined with T d = 18 K would still result in cold gas, T K = 4K. A distinct AT 
for these points is an unlikely explanation for their high r 240 ratio. 

The optical depth of the 13 CO J = 1 — ► line is high for these points, while T K is 
still well above 3 K. This allows r 240 to be high, while obviating the need for cold gas. 
Given that the models require component 1 to be cold for these points (i.e., T K = 3 
to 10 K) and component to be 18 K, then the models already require at least some 
of this gas to be optically thick in the 13 CO J = 1 — > line and some to be optically 
thin. So, if the observations imply a high 13 CO J = 1 — > optical depth for all of the 
gas, then the models would be in error. The gas density is near or above the critical 
density of the 13 CO J = 1 — > transition, so the line is close to LTE. Consequently, 
the 13 CO J = 1 — > 0/ 12 CO J = 1 — ► intensity ratio is a good estimate of the optical 
depth of 13 CO J = 1 — > 0. This ratio for these points is about 0.3 — significantly less 
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than unity. Therefore, at least some gas is indeed optically thin in 13 CO J = 1 — > 0, in 
agreement with the models. 

A lower 13 CO abundance or a higher dust-to-gas ratio for these points could account for 
the high r 240 values. Accordingly, X( 13 CO) would be factors of 2 to 3 lower or x d would 
be similar factors higher. There is no observational evidence for large ch anges in the 



13 CO abundance or the dust-to-gas mass ratio. Nevertheless, the models of IWarin et al. 



240 

n 



( 119961 ) suggest that selective photodissociation of 13 CO reduces its abundance by factors 
of 2 or 3 near molecular cloud surfaces. If the triangular cluster of points in the r 
versus T dc plot were affected by this reduced X( 13 CO), then correcting to the "normal 
abundance would bring these points down by factors of 2 to 3. This would increase the 
lower temperature limit of the cold gas. The bend in the model curve from vertical to 
horizontal occurs more or less at the bottom of the triangular cluster and more or less 
for a component-1 temperature of 7K. Consequently, correcting for a possible reduced 
13 CO abundance could possibly bring the lower temperature limit up to slightly less 
than 7K. 

A similar possibility would be an emissivity enhancement of the dust: the dust mass 
absorption coefficient, i.e. k v at 240 /xm, would be unusually large, by factors of 2 
to 3, for these points. For example, a resonance in the dust absorption spectrum 
at 240 /mi could produce an increase in k v at 240 /im. However, an increase by the 
same factor would be needed at 140 /im in order to maintain Td c at about 18 K. This 
resonance feature would be at least 100 /im wide. This is very unlikely. Instead of 
a resonance, dust with o pacitie s 2 to 3 times higher than normal dust at all FIR 



wavelengths is sufficient. iDwekl (120041 ). for example, discusses dust grains with far- 
IR opacities orders of magnitudes higher than the classical silicate, graphite grains. 
Mixing a very small portion of such grains with classical grains could give a mix with 
an effective far-IR opacity easily factors of 2 to 3 higher than considered here. However, 
the observations require the opacity to increase as the r 240 value increases towards the 
top of the triangular cluster. Increasing opacity at far-IR wavelengths increases cooling 
as well. Consequently, there would be an obvious overall trend to lower and lower 
140 /im/240 /iin color temperature as r 240 increased — resulting in a triangular cluster 
whose peak would be noticeably skewed towards lower colour temperatures. This 
is not observed. Increasing the far-IR emissivity while keeping the 140/zm/240//m 
color temperature constant requires having another warmer component mixed in that 
dominates the dust continuum emission. But this is nothing more than the original 
two-component models that have been used up to this point. Also, grains with such 
high far-IR emissivities are more likely to be responsible for such low dust temperatures 
(i.e. about 3 - 5K) rather than rule them out. 
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• Yet another possibility is that points towards the peak of the triangular cluster are 
dominated by a component with a large value of r 240 and Td c around 18 K. In terms of 
the two-component models discussed in the current work, this is equivalent to increasing 
the parameter cq for the peak of the triangular cluster. 

It seems that at least some alternatives to cold dust and gas may exist. Overall, however, 
the evidence is far from convincing. In addition, simply abandoning the basic assumption 
leads to ad hoc interpretations. This is probably unjustified given the success of the basic 
assumption, and its associated models, at accounting for the overall trend in the data; the 
basic assumption and the models should not be so lightly discarded. 

Can we explain the data without such cold dust and gas and still use the models and 
the basic assumption? In other words, is there some unexplored region of parameter space 
that permits a higher lower limit on the dust and gas temperature of component 1? The 
answer is yes, but the changes are not dramatic. The lower limit to T K1 and T^i depends on 
the position of the "component-0 point" . That is, the physical parameters of component 
are constant in every respect and therefore represent a single point in the r 240 versus Ta c 
plot. In contrast, the gas and dust temperatures of component 1 vary spatially, while the 
other parameters are held constant, and therefore the parameters of component 1 represent a 
locus of points in this plot. The position of the component-0 point is just above the vertical 
section of the plotted two-component model curve. When the r 240 of the component-0 point 
is higher than the apex of the triangular cluster, then the lower limit on Tdi is higher. In 
fact, this lower limit is roughly given by Tdi on the model curve at the apex position. Hence, 
the relevant parameter space area is where the component-0 point has higher r 240 . This is 
easily accomplished by reducing the 13 CO abundance of component only by a factor of 
2. Doing this, and refitting the two-component models, raises the component-0 point from 
r 240 ~ 65 to about 110 MJy • sr _1 • (K • km • s -1 ) -1 . However, simply keeping the 13 CO 
abundance fixed at its usual value and restricting c to be > 1 results in even higher r 240 for 
the component-0 point. In spite of the more promising position of this point, T<ji at the apex 
of the triangular cluster changes only by about 0.5 K. If we truncate the model curve above 
the Tdi = 5K point (i.e. remove the temperatures lower than this), and keeping AT > OK 
to keep T K1 > 5 K as well, then the xl increases from 5.6 to 8.1, an increase excluded by the 
F-test at a confidence level of more than 90%. Accordingly, raising the component-0 point 
to higher r 240 merely stretches the vertical section of the curve between Tdi — 2.8 K and 
Tdi — 5 K. Consequently, we only increase the lower limit on Tdi from about 3 K to about 
3.5 K. In short, there is no compelling evidence that rules out cold dust and gas. 



Previous evidence for cold dust or gas is not compelling eithe r , but does nonethe- 
less come from a wide variety of observations (e.g., iReach et al.lll995l ; iMerluzzi et al.lll994j ; 
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Ristorcelli et al.lll998l). In addition, there is evidence for gas that had be en previously unde 



tected fe.g-. lReach et alll998l : ICuillandre et al.ll200lh . iReach et al.lfll995f ) used COBE/FIRAS 



continuum data with observed wavelengths from 104 /im to 2 mm to infer a widespread cold 
component with dust temperatures 4 to 7K. This c omponent is found at all Galactic lati- 
tudes from the Galactic plane to the Galactic poles. Lagache et al. ( 19981 ) re-examined the 



FIRAS data and concluded that the cold componen t oflReach et al.l ( 19951) was not needed 



the coldest component necessary was at about 15 K. iFinkbeiner et al.l (119991 ). in yet another 
examination of the FIRAS data, co ncluded that the colder component has a temperature 
of about 9 K. iMerluzzi et al.l (119941 ) found a cold component with a temperature of either 
15 K o r 7K, depending on w hether the spectral emissivity index, /3, was 1.1 or 2, respec- 
tively. iRistorcelli et al.l (119981 ) armed with continuum observations in four wavelength bands 
(i.e., 180-240, 240-340, 340-560, 560-1050 /im) discovered a "cold condensation" close to the 
Orion Nebula with a temperature of 12.5±3 K. This is not as cold as the 3 to 10 K material 
discussed here, but nevertheless shows that dust with temperatures significantly lower than 
the 18 K expected for dust heated primarily by the general ISRF is possible. In contrast to 
these previous papers, the current paper infers a cold component (i.e. T~ 3 to 10 K) without 
benefit of long-wavelength (i.e. 1 mm) continuum data. The long-wavelength data used 
here is the 2.7mm 13 CO J = 1 — > spectral line in emission. 

If we accept for the moment that a cold component with temperatures of 3 to 10 K does 
indeed exist within the Orion clouds, then the obvious question is how can such cold gas and 



dust e xist without being strongly affected by the general ISRF or local stars? IReach et al. 



( 119951 ) discuss a number of possible explanations in the context of cold dust throughout the 
Galaxy, especially in high-latitude clouds. We revisit some of the proposed explanations of 
Reach et al.l (119951 ). but in the context of the Orion clouds: 



Shielding from the Interstellar Radiation Field. As lReach et al.l (119951 ) point out, atten- 
uating the he ating rate by a facto r of 10 3 requires a minimum absorption equivalent to 
A v = 20 mag (IMathis et al.lll983l ). Given that the radiation field is proportional to 
and that T d = 3 to 10 K is the cold dust temperature range, then attenuation factors of 
roughly 30 to 5x 10 4 are necessary. This then requires mini mum absorptions eq uivalent 
to A v ~ few mag to A v considerably more than 20 mag (IMathis et al.lll983l ). These 
correspond to column densities of a few x 10 21 to more than 3 x 10 22 H nuclei- cm" 2 . 
Over much of their area, the Orion clouds have N(H) closer to the former value. Con- 
sequently, shielding might explain the Tj = 10 K material, but is unlikely to explain 
the really cold material (Td ~ 5K). Also, since most of this cold material is on the 
cloud edges, shielding probably will not account for the cold dust and gas suggested 
by the models. 
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Fractal Grains. As previously mentioned, grains with enhanced sub millimete r /FIR 
emissi on relative to UV/vis ible absorption, li ke the iron needles discussed by iDwek 
( 120041 ) or fractal grains (see lReach et al.lll995l . and reference s there in), can have tem- 
peratures near that of the cosmic background. iReach et al.l (119951 ) state that fractal 
grains may have greatly reduced mass for the given absorption cross-section; the 60% 
upward correction to the single-component mass would be revised considerably down- 
ward. In addition, fractal grains may also have a reduced volume for the given geomet- 
ric cross-section, permitting AT ~ 0, as discussed previously (see also Appendix IA~|) . 

Very Large Grains. Dust grains larger than the FIR wavelengths t hat t hey emit do so 
very efficiently and, therefore, cool very efficiently. IReach et al.l (119951 ) find that the 
size dis tribution of these gr ains must steepen beyond that of the power-law index of 
—3.5 of iMathis et al.l (119771 ) to explain their observations. This steepening prevents a 
dust-to-gas mass ratio much higher than observed; the power-law index of —3.5 out 
to a maximum grain radius of 100 fim increases the dust-to-gas mass ratio by more 
than an order of magnitude. In addition, the current work requires the observed grains 
to have a large geometric cross-section to volume ratio, thereby ruling out such large 
grains being the bulk of the cold dust. 



Long-Wavelength Emissivity Enhancement. IReach et al.l (119951 ) suggest that an en- 
hancement of the continuum emissivity at A ~ 800 fim could explain their observations 
instead of some hypothetical cold dust. Obviously, this explanation does not apply to 
the current paper because the models have been applied to 240 fim continuum data 
and 2.7 mm spectral line data. If the cold material pre dicted by the mod els for the 
Orion clouds is indeed ubiquitous in the Galaxy, then the IReach et al.l (119951 ) proposed 
enhancement of the long- wavelength continuum emissivity cannot rule it out. 

Of the explanations given above, fractal grains may be the most feasible. 

Thus the observations neither completely exclude nor strongly support the existence of 
cold dust and gas (i.e. T~ 3 to 10 K) in the Orion clouds nor a widespread presence in the 
Galaxy as a whole. The least unlikely alternative to dust and gas at temperatures of 3-5 K is 
probably that mentioned in Subsection 13.21 the additional emission of dust grains associated 
with atomic hydrogen. Nevertheless, it seems likely that cold material with temperatures of 
7 to 10 K does exist in the Orion clouds. If cold dust and gas exist in the ISM in general 
with temperatures of ~ 3 to 10 K, then fractal dust grains or iron needles may be the most 
credible reason and the 60% upward correction to the Galactic ISM's mass due to this cold 
material might be revised downward substantially. This downward revision is also necessary 
if the material at temperatures of 3 to 7 K is really an artefact due to the other effects just 
described. 
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4.6. The Millimeter Continuum to 13 CO J = 1 —> Ratio as a Temperature 

Diagnostic 

Given that we can now characterize the 1^(240 /mi) /I ( 13 CO J = 1 — > 0) ratio as a 
function of physical parameters like the temperature, continuum and 13 CO J = 1 — > ob- 
servations can constrain the temperature of the dust and molecular gas. In particular, using 
only ground-based observations to achieve such a constraint is advantageous. Consequently, 
we examine one particular representative case: using the 1(1300 /im)/I( 13 CO J = 1 — > 0) 
ratio — hereafter, r 1300 — to estimate the dust temperature, T,j. 

To this end, the simulations described previously in Paper II were used to generate 
a model 1300 jum continuum map. The adopted mass absorption coefficient at 1300 /im, 
K„(1300/xm), was determined from the adopted ^(lOO/xm) of 40 cm 2 • g _1 and scaling by 
v p for (3 = 2, yielding k„(1300 /an) = 40 cm 2 ■ g" 1 x (100/1300) 2 or 0.24cm 2 • g~ x . To adopt 
a reasonable noise level, we note that current bolometers can achieve an rms noise value of 
2 MJy-sr" 1 for a 1-hour integration in a 15" beam at wavelength of 450 /im (D. Hughes, priv. 
comm.). Assuming a noise proportional to A -2 and, depending on the integration time and 
spatial averaging of the observations, a noise level of 0.075 MJy • sr _1 at 1300 /im is possible. 
Therefore, this rms noise level is adopted for the simulated 1300 /im map. The simulated 
1300 /im continuum map was divided by the simulated 13 CO J = 1 — > map to produce the 
r 1300 map. Figure H] shows the resultant simulated r 1300 values plotted against the simulated 
Tdc values — the 140 /im/240 /im color temperature. This color temperature is, of course, 
not directly relevant to 1300 /im observations, but is included in Figure H] to permit easy 
comparison with the standard r 240 versus T dc plots used throughout the current paper. 

The upper panel of Figure 0] shows that r 1300 is not useful as a temperature diagnostic 
for dust color temperatures of ~ 15 to 30 K; given that the vertical spread of the data points 
is about the size of the vertical error bars, we cannot unambiguously associate each r 1300 
with a single T dc . The temperature dependence is weak in this range, because both the 
continuum and line observations are close to the Rayleigh- Jeans limit. This implies that 
neither component's emission will overwhelmingly dominate over the other's for a larger 
range of temperature differences between the two components. In contrast, at A = 240 /im 
and in the ~3 to 30 K range of dust physical temperatures, component 1 goes from being 
overwhelmingly dominated by to overwhelmingly dominating the component-0 emission. At 
A = 1.3 mm, component- 1 goes from being dominated by only factors of a few to dominating 
by only a factor of ~2. Accordingly, combined emission of the two components has a narrow 
range of brightnesses. 

In contrast, the lower panel of Figure @] demonstates that r 1300 can indeed be a tem- 
perature diagnostic. This plot shows the model curve up to temperatures of 200 K, and the 
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simulated map data points in the lower left corner. The curve goes from nearly flat (i.e. 
slope near 0) for Td c ~ 15 to 30 K to a power law with 

r i300 00 ^d' 4 (3) 

above a threshold of 

r^^O.SMJy-sr" 1 • (K-km-s -1 ) -1 and T d >50K . (4) 

(Note that Tj can serve in place of Td c for high temperatures, because the effect of com- 
ponent is negligible in this limit.) Therefore, r 1300 is a useful temperature diagnostic for 
temperatures above about 50 K or r 1300 above ~ 0.5 MJy • sr" 1 • (K • km • s -1 ) -1 . Also, above 
this threshold the temperature sensitivity increases with increasing temperature. This is not 
quite as strong as the T 2 rise expected in the LTE, high-temperature limit, but is definitely 
much better high-temperature sensitivity than for the ratio of two rotational lines of CO or 
for the ratio of two continuum bands. 

Naturally, the real situation is considerably more complicated than simply reading the 
temperature from this one curve that may apply only to the Orion clouds. As we recall, there 
are a number of physical parameters that affect this curve: Tdo, ^r, Co, n c0 , n c i, AT, n u 
(or (3), x d , and N(HI)/N(H 2 ). Before we examine how variations of these parameters effect 
the model curve, the curve in Figure H] is adopted as the nominal curve. The Td inferred 
from an observed r 1300 using this nominal curve we call the nominal T<j. Then the question is 
how the true Td differs from the nominal Td due to physical parameters differing from their 
nominal values (see the "Input Values" column of Table 1 of Paper II). A reasonable goal is a 
temperature estimate within a factor of 2 of the true temperature for those "hot spots" where 
r 1300 is above the threshold value of 0.5 MJy-sr -1 • (K-km-s -1 ) -1 . If we vary the component-0 
parameter values (i.e. the first four parameters listed above), then the true Td stays within 
a factor of 2 of the nominal Td in most cases. Specifically, varying Tdo from values lower 
than nominal to double this, Co^- from factors 10 or more lower (since this represents the 
optically thin limit), and n c o from 10cm -3 to 10 5 cm -3 and higher (since this represents the 
LTE limit) changes r 1300 by only 5-7%. Equivalently, the true Td differs from the nominal 
value by 7-10%, well within the desired factor of 2. There are two potential difficulties. One 
is if co|^f- is two orders of magnitude higher than the nominal value, then the true r 1300 is 80% 
lower than its nominal value. In practice this would not be a problem because the observed 
r 1300 would be low enough that this observed position would not be identified as a hot spot; 
there would be no false positives in the search for hot spots. Nevertheless, any false hot spot 
is easily identifiable with supplemental observations of the J = 1 — > line of C 18 or 12 CO; 
the C 18 0/ 13 CO intensity ratio would identify the position as having high 13 CO J = 1 — > 
opacity. The 13 CO/ 12 CO intensity ratio could be used similarly, but C 18 0/ 13 CO would be 
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more reliable given that the corresponding abundance ratio is within an order of magnitude 
of unity. Another potential difficulty is that, even though varying n cQ over many orders of 
magnitude has a negligible effect on r 1300 for T dc > 50 K, n c0 as low as 10 cm -3 can raise the 
spur of data points at T dc ~ 18 K (see Figure Hj) by about 70%. Accordingly, the threshold 
listed above, expressions (jl]), is raised to 

r 1300 >0.85MJy • sr" 1 • (K-km-s -1 ) -1 and T d >60K . (5) 

Densities as low as 10 cm -3 are probably unlikely, so threshold (j5J) is probably unnecessary. 
Nevertheless, it can be used if extra caution is desired. 

Varying the component- 1 parameters and n cl results in the true T d being up to 
a factor of 5 smaller than the nominal T d . The value was increased by an order of 
magnitude or decreased similarly (or more because this is the optically thin limit). At the 
same time, the n c \ value was decreased to 100 cm -3 (again assuming that densities as low 
as 10 cm -3 are unlikely) and increased to 10 5 cm -3 (or more because this is the LTE limit). 
Over most of this parameter space the true T d stayed within a factor of 2 of the nominal 
Td. However, the combination of n& as low as 100 cm -3 and |^ an order of magnitude 
larger gives a true Td up to a factor of 5 smaller than the nominal Td. This low temperature 
case can be identified by observing the 13 CO J = 2 — ► 1 line. In the nominal case, the 
ratio I( 13 CO J = 2 -> l)/I( 13 CO J = 1 -> 0) varies between 1.9 and 3.4 for T K = 50 to 
200 K. In this high-§£, low-n cl case, I( 13 CO J = 2 -> l)/I( 13 CO J = 1 -> 0) is 0.5 to 0.7. 
Supplemental continuum observations help as well. The L,(450 /j,m)/l v (1300 /xm) ratio, for 
example, is sensitive to temperatures up to T d — 60 K for a flux ratio uncertainty of 20%. 
Consequently, the observed 1^(450 /xm)/I I/ (1300/xm) ratio places a lower limit on T d . 

Considering variations in AT, the determination of Td is remarkably insensitive to such 
variations. For example, AT could range from —80 K to +80 K and the true Td is still within 
about 20% of the nominal Td. Of course, the T K would be quite different from Td in that 
case. For AT = —80 K, the T K is even higher than a nominal Td that was already high; this 
T K is still within a factor of roughly 2 of the nominal Td- For AT = +80 K, T K is quite small. 
The I( 13 CO J = 2 -> l)/I( 13 CO J = 1 -> 0) ratio is low in this case, but it is difficult to 
distinguish this case from the low-Tj, low-density case described in the previous paragraph 
using only this CO-line ratio. Using the L,(450 /zm)/I„(1300/zm) ratio in conjunction with 
the I( 13 CO J = 2 — > l)/I( 13 CO J = 1 — > 0) may provide sufficient information to distinguish 
this high-T d , low-T K case from the low-Td, low-density case. 

Dealing with anomalous values of the remaining parameters - - k v (or 0), x d , and 
N(HI)/N(H 2 ) — is problematic but still possible. If a source (a molecular cloud or part 
of another galaxy) is mapped in L,(1300 /im) and 13 CO J = 1 — > 0, then the majority of the 
positions probably have roughly constant values for these parameters; the majority of these 
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positions would also have similar values of r 1300 . For example, it is easy to imagine that the 
molecular clouds in another galaxy have an unusually high or low dust-to-gas ratio, i.e. high 
or low x d , compared with such clouds in our Galaxy. Or we can imagine a cloud that has a 
non-negligible layer of atomic gas (i.e. non-negligible N(HI)/N(H 2 )). In cases such as these, 
the r 1300 value itself would not be important, but that value in relation to the "average" r 1300 
values for the observed source. The hot spots are identified as those with r 1300 values large 
compared to the typical value for the source. If, instead, such hot spots are not due to ele- 
vated temperatures but unusual k V) x d) and/or N(HI)/N(H 2 ) values, then such "hot spots" 
would still be interesting: they represent positions with unusual properties. Specifically 
identifying those unusual properties would involve the supplemental observations described 
above — using the I( 13 CO J = 2 -> l)/I( 13 CO J = 1 -> 0) and 1,(450 /W) /I„( 1300 fjm) 
ratios — and other observations when possible. 

Addressing the various points mentioned above, the following is a plausible observing 
plan: 

1. First map the 1300 /im continuum and the 13 CO J = 1 — > line. The majority of 
points establish the typical r 1300 value for the source. Any points with high r 1300 values, 
above some threshold that applies to the source, would be hot spots to be ear-marked 
for further observations and study. The r 1300 values for these hot spots will provide a 
temperature estimate to within a factor of 2 in the majority of cases. 

2. Map the 13 CO J = 2 -> 1 line. The I( 13 CO J = 2 -> l)/I( 13 CO J = 1 -> 0) ratio can 
help confirm whether the hot spots are indeed hot. 

3. Map the 450 /im (or shorter) continuum. The 1^(450 jum) /I„(1300 /im) ratio is a further 
check on the temperature and can, in some cases, check whether AT = or not. 

4. Map the C 18 J = 1 -> line. The I(C 18 J = 1 -> 0)/I( 13 CO J = 1 -> 0) ratio can 
confirm that 13 CO J = 1 — > is indeed optically thin, thereby further checking the 
temperature inferred from the r 1300 of the hot spots. Even though the 13 CO J = 1 — > is 
expected to be optically thin on multi-parsec scales, as in the Orion clouds, verification 
of this can rule out other possibilities. 

All of these observations are possible from the ground. 

Obviously, this is only one proposed observing plan out of many possibilities. Similar 
plans could be devised using other continuum wavelengths and other rotational lines. Two 
such alternatives would be the L,(2700 /zm)/I( 13 CO J = 1 — > 0) ratio or the 
L,(1300 /xm)/I( 13 CO J = 2 — > 1) ratio. This is re-discovering the line-to-continuum (simi- 
larly, the continuum-to-line) ratio or the equivalent width, but for millimeter- or 
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submillimeter-wave molecular lines. Such equivalent widths are commonly used at visible, 
i nfrared, and cen t imeter wavel engths for inferring the physical properties of HII regions 
( 10sterbrocklll989l ; ISpitzerl Il978l ) . Likewise, the equivalent widths of millimeter- wave molec- 
ular lines can provide important physical insights into molecular clouds (i.e. H 2 regions). 

Nevertheless, the 13 CO J = 1 — > line has a big advantage over higher rotational lines 
of 13 CO: relative insensitivity to physical parameters like density and column density per 
velocity interval. This is why the r 1300 ratio can often predict to within a factor of 2, 
despite varying n(H 2 ) and N( 13 CO)/Av by orders of magnitude. At the very least, this 
method potentially places a realistic and interesting upper limit on dust and molecular gas 
temperatures. 



5. Conclusions 



Far-infrared continuum data from the DIRBE instrument aboard the COBE spacecraft 
were combined with 13 CO J = 1 — > spectral line data from the Nagoya 4-m telescope 
to infer the large-scale (i.e. ~ 5 to ~ 100 pc) physical conditions in the Orion molec- 
ular clouds. The 140 /im/240/zm dust color temperatures, Td c , were compared with the 
240/im/ 13 CO J = 1 — > intensity ratios, r 240 , to constrain dust and molecular gas physical 
conditions. In addition, such a comparison provides valuable insights into how the ratio of 
FIR/submillimeter /millimeter continuum to that of a 13 CO (or C 18 0) rotational line can 
constrain temperature estimates of the dust and molecular gas. For example, ratios of rota- 
tional lines or ratios of continuum emission in different wavelength bands often cannot place 
realistic upper limits on gas or dust temperature, whereas the continuum-to-line ratio can 
place such limits. 

Two-component models fit the Orion data best. One component has a fixed-temperature 
and represents the gas and dust towards the surface of the clouds and is heated primarily 
by a kiloparsec-scale interstellar radiation field, referred to here as the general ISRF. The 
other component has a spatially varying temperature and represents gas and dust towards 
the interior of the clouds that can be both shielded from the general ISRF and heated by 
local stars. The model results and their implications are as follows: 



1) The inferred physical conditions are consistent with those deri yed from the lar g e-scal e 
observations of the J = 2 — > 1 and J = 1 — > lines of 12 CO by lSakamoto et al.l (119941 ) . 



At least two gas (dust) temperatures are needed on the majority of sightlines through 
molecul ar clouds for reliabl y estimating column densities. This is supported by the 
work of ISchnee et al.l (120061 ) . 
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3) The dust-gas temperature difference, Td — T K or AT, is K to within 1 or 2 K. If this 
result applies more generally to the Galactic-scale molecular ISM, except for unusual 
regions such as the Galactic Center, then there are a number of implications: 

— Dust-gas thermal coupling is factors of 5 to 10 stronger than has been previously 
assumed. Such factors may be due to the distribution of dust grain sizes and 
grains with larger cross-section to volume ratios than that of a simple sphere. 

— Galactic-scale molecular gas temperatures are closer to 20 K than to 10 K, because 
the emission from the CO rotational lines, even the optically thick 12 CO J = 1 — > 
line, does not fill the beam within the velocity interval about the line peak. 

— This CO emission that does not fill the beam provides a better explanation of 
the NfhhjVlf CO) conversi on factor or X-factor. Discussion of this is deferred to 
a later paper JWallll2007bh . 

— Having AT nearly constrains which mechanisms explain dust grain alignment 

in the ISM. A negligible dust-gas temperature difference rules out t he Davis- 

Greenstein alignment me chanism, but not other possible mechanisms (see lLazarian et al 



19971 ; I Abbas et al.l 12004 . and references therein) . 



4) Roughly 40-50% of the ISM in Orion is cold (i.e. 10 K) to very cold (i.e. down to 
3 K) dust and gas. Accordingly, there is roughly 60% more gas and dust in Orion than 
inferred from simple one-component models. This may also imply a similar increase in 



the estimated mass of entire Galactic I SM. Fractal dust grains (see lReach et al.lll995 



and references therein) or iron needles (jDwek!l2004l ) may explain the low temperatures 
of this gas and dust and, at the same time, may account for the high dust-gas thermal 
coupling needed to explain AT ~ OK. Nevertheless, alternative explanations that do 
not require cold dust and gas cannot be ruled out; the least unlikely of these other 
explanations is a contribution to the 240 fim continuum emission of the dust associated 
with atomic hydrogen. The data suggest that the effect of the H I-associated dust is 
negligible, but still might permit raising the lower temperature limit of this cold gas 
and dust from 3 to 5 K. 



The model parameter values derived from the fits to the r 240 versus T^c plots were 
used to create simulated 1300 /im continuum and 13 CO J = 1 — > line maps. These 
simulated maps tested whether the millimeter continuum to 13 CO J = 1 — > line inten- 
sity ratio could constrain temperature estimates of the dust and molecular gas. The ratio 
L,(1300 /zm)/I( 13 CO J = 1 — ► 0), or r 1300 , was found to estimate the dust temperature to 
within a factor of 2 in most cases, provided that r 1300 was higher than a threshold level of 
0.5 MJy ■ sr _1 ■ (K ■ km ■ s -1 ) -1 . Supplemental observations of the 13 CO J = 2 — > 1 line and 
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shorter wavelength continuum would confirm the high temperatures in these high-r 1300 "hot 
spots" . The results here can be easily generalized to other continuum wavelengths and other 
rotational lines, even permitting interpretation of millimeter and submillimeter molecular 
line equivalent widths. And this is entirely possible with only ground-based observations. 

The full potential of using millimeter continuum and 13 CO (or C 18 0) rotational line 
comparisons has yet to be realized. 
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A. The Effect of Grain Size on Gas-Grain Thermal Coupling 



Burke fc Hollenbachl (119831 ) describe the heat transfer rate from the gas to the dust (or 



vice versa as 



A„ 



n„„n a 

H2 gr gr 



8kT T . 



7THl t 



a T (2k)(T K -T d ) 



(Al) 



where n H2 and n gr are the number densities of hydrogen molecules and dust grains, respec- 
tively, a r is the grain geometric cross-sectional the gas kinetic temperature, T^ is 
the dust grain temperature, m H2 is the mass of the hydrogen molecule, and a T is the accom- 
modation coefficient. This coefficient is a measure of how well the temperature of the gas 
particles that have collided with the dust grains accommodate to that of the grain surface 
( Burke &: Hollenbach 19831). (Note that the n and m H2 in this expression replace the n H 
and m H in expression 9 of iBurke fc Hollenbach ( 1983 ) because the colliders considered here 
are hydrogen molecules.) If we simplistically assume identical grains, then 



or 



n . 



n m 

gr gr 



• E d I1 H2 m H2 

m„ 



(A2) 
(A3) 



-39- 



where m gr is the mass of a single dust grain. Assuming spherical grains of uniform density, 
p yields 



and 



m„ 



a 



3™ Pgr 



(A4) 



(A5) 



where a is the grain radius. Expression ( 1A4I) substituted into ( 1A3I) gives 



n 



and multiplying the above by ( 1A5I) yields 



£ X d n H2 m H2 

4 vra 3 p gr 



3 X d n H2 m H2 



(A6) 



(AT) 



Note that expression (IA7I) is nearly identical to expression (5) of iGoldsmithl (120011 ); he also 
included a Q-correction factor that is not necessary here because we are only concerned 
with the geometric cross-section and not the absorption cross-section. Note also that ( 1A7I) 
is proportional to the grain cross-section to volume ratio. (In reality, it is the grain cross- 
section to mass ratio that is relevant, but, for a uniform grain density, this is equivalent to 
a proportionality to the cross-section to volume ratio.) Substituting (1A7|) into (lAlj) yields 



A 



gd 



3 
2 



8A; 3 



m T 



7T 



Oi„ 



X,, 



ap K 



n 



(T K - T d ) t; 



(A8) 



Putting in the physical constants in cgs units gives 

x 



A 



gd 



7.11 x 10 



" 36 a, 



ap 



' n 2 

H2 



(T K - T d ) T^ 



(A9) 



Goldsmith! (120011 ) adopted the following values for the dust pa r amete rs 



x , 



2g • cm 3 , and a = 
a = 1.7 x 10~ 7 cm" 



1.7 x 10~ 5 cm- 3 . 



0.01, Ps 



[Note that IGoldsmithl (120011 ) apparently adopted 
" 3 , but this is near the lower limit of the dust grain size range (see 
Desert et al.lll990l ). Also, and more importantly, that value of a is inconsistent with the 
numerical coefficients in the expressions for A g d. Finally, a = 1.7 x 10~ 5 cm~ 3 was the 
actual intended dust grain radius (Goldsm ith, priv. comm.).} T ne re commended value of 
the accommodation coefficient is a T = 0.3 (IBurke fc Hollenbachlll983l ). Consequently, 



A„ d = 2.0 x 10 
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,(T K -T d ) 



10 K 



(AlO) 



as per 



Goldsmith! (1200 if ). 
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However, as stated in Section l4~4t A g a must be factors of 5 to 10 larger to explain the 
observations. This is achieveable using a more realistic treatment of grain sizes in deriving 
A gd . Specifically, a range of gr ain sizes must be con sidered instead of simplistically adopting 
a single radius. For example, iDesert et al.l (119901 ) suggest that a = 15 to 110 nm for the 
big thermal equilibrium grains. Consequently, the 170 nm adopted by iGoldsmithl (120011 ) is 
clearly too large; instead, some appropriately weighted mean of 1 5 nm and 110 nm is the 
most realistic choice for the a in expression (1A9j) . From the work of iMathis et al.l (119771 ) we 
know that 



n (a) 



A/q a 



-3.5 



(All) 



where n (a)do is the number density of grains with radii between a and a + da and k is a 
normalizing constant. The k is found from 

I max 

a' 3 - 5 da , (A12) 



n (tot) 



with n (tot) as the number density of all the grains and corresponds to the n gr in (1A2j) and 
i3~l) . From (1A12[) . k can be expressed as 



k n 



n_(tot) k x 



(A13) 



k, is another constant of the distribution and depends on a and a . Its exact form is 

1 Sr mtn max 

not relevant to the derivation here, but is nonetheless given below for completeness: 

(a a ) 2 - 5 , s 

2 V mm max/ f A, 14:1 



k, 



a 2 - 5 — a 2 - 5 



To include the effect of the range of grain sizes on A gd , the n gr cx gr in equation (I A II) must 
be replaced with n gr (tot)(<j gr (a)) a , where (cr gr (a)) a is the grain cross-section after averaging 
over the size distribution. Thus, 



n g r( tot )( ff g ») ( 



n gr( Q ) °"gr( a ) da 



,0.5 



,0.5 



2nk 1 n (tot) 



(a . a ) - 5 

V mm max / 



(A15) 
(A16) 



where equations (IA11I) . (IA5I) . and (IA13j) were used. The n (tot) on the right side must now 
be expressed in terms of the molecular gas density, n H2 . This can be done using an expression 
analogous to that of ( 1A3I) that uses the range of dust sizes: 



n gr (tot)(m gr (a)). 



a 'd ri H2 m H2 



(A17) 



or 
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/■a 

/ max 

/ n gr (a)m gr (a)da = x d n H2 m H2 



Applying equations (jA~4j) . flAlip . and (IA13[) yields 

y/'/^NU a~ a5 da = x d n H2 m 

Integrating and solving for n (tot) gives us 



MtOt) = 3X d n H2 m H 2 (Alg) 

srV ; 8vrp jfe. a - 5 -a - 5 ^ 

' g r 1 v max r, 



rn i n ' 



Substituting ( 1A18j) into the right side of ( 1A16I) yields 

n» t (tot)<^(«)>a = I^T 1 > ( Al9 ) 

where 

a = (a . a ) a5 . (A20) 

e/j V mtn max J \ I 

Equation ( 1A19I) replaces the n gr cr gr that appears in ( 1A1I) . yielding an expression nearly iden- 
tical to ( 1A8I) and (1A9jl . except that a is replaced with a eff . Therefore, the relevant grain 
radius in the expressions for A g( j (e.g., equation IA9I) is the geometric mean of the minimum 
and maximum grain sizes (i.e., equation IA20j) . 

The a and a should be those for the big grains, rather than for the full range of 

mtn max O O 7 O 

dust sizes that also inclu de the VSGs (very small grains) and the PAHs (polycyclic aromatic 



hydrocarbons) (e.g., see lDesert et al.lll990l ). The big grains are in thermal equilibrium and 



are the grains obse rved with the 140 /im and 240 fim DIRBE observations. According to 



Desert et al.l (119901 ). the a min and a max are 15 and 110 nm, respectively, for the big grains. 
This gives a eff = 41 nm and increases A g d by a factor of 4.2. The full range of sizes over 
all grains, i.e. a = 0.4 nm and a = 110 nm, results in a ,, = 6.6 nm and A„h is 

O i min max 7 ej j B u 

increased by a factor of 26. However, there are at least two problems with using the full size 
range of grains. One is that, as mentioned above, only the big grains are relevant to the 
observations discussed here. The second is that the treatment above implicitly assumes that, 
within the grain size distribution, only the grain size changes; other grain properties, such 
as grain density and shape, are assumed constant despite changing grain size. For example, 
going from (1A17I) to (1A18|) assumes that p gr is independent of a. This is likely to be a bad 
ap proximation for the full size range, especially when grain type varies with grain size (e.g., 



sec 



Desert et al.lll990l ). Therefore, the factor of 4.2 increase in A gc j is appropriate when only 



considering the big grains. 

A few corrections should be considered before using that factor of 4.2. Given that 
only the big grains were used, we must correct for not using the full population of dust 
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grains. Specifically, the dust-to-gas mass ratio used must be replaced by th e mass ratio 



of the dust in big grains to that of the gas. The grain densities and sizes in IDesert et al. 



( 11990l ) imply that the big grains represent 76% of the mass of the dust. However, a more 
appropriate accommodation factor is necessary for the low dust and gas temperatures tha t 
considered here. This suggests that a T is 0.4 instead of 0.3 (see lBurke fc Hollenbachlll983l ). 
Accordingly, increasing the accommodation factor while decreasing the dust-to-mass ratio by 
similar amounts gives an overall correction of nearly unity (in fact it is about 0.9). Another 
possible co r rectio n, or at least uncertain ty, is the density ad opted for the big grains. While 



Goldsmith! (120011 ) adopts p = 2 g-cm 3 , IDesert et al.l (Il990l ) use 3 g - cm 3 for the big grains 



This latter density brings the 4.2 factor down to about 3. 

In conclusion, considering a realistic range of grain sizes increases the gas-grain thermal 
coupling by factors of about 3 to 4. 
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Fig. 1. — Plots of the 140 /xm/240 //m dust color temperature versus the ratio of the atomic 
gas to molecular gas column densities, N(HI)/2N(H 2 ), appear above. The upper two panels 
show these plots for the one-component, LVG models. The lower two panels are for the 
two- component, two-subsample, LVG models. The panels on the left include the error bars, 
while the panels on the right exclude the error bars. The curves in the panels on the right 
represent hypothetical cases where the dust associated with the molecular gas has one fixed 
temperature for all lines of sight and the dust associated with the atomic gas has another 
fixed temperature. The lower curve in each of the right panels assumes that the dust in 
the molecular gas has T d = 16.5 K and the dust in the atomic gas has T d = 22.5 K. The 
upper curve in each of the right panels assumes T d = 27 K and 22.5 K for the dust associated 
with molecular and atomic gas, respectively. The plots only include those pixels with the 
intensities above the 5-cr level in I„(140 /xm), I„(240 /iin), I( 13 CO) and above the 3-a level in 
I(HI) simultaneously. 
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Fig. 2. — The is the equivalent of Figure [U but for the simulated data. The HI layer has 
a constant column density of 5xl0 20 if atoms ■ cm -2 and a constant dust temperature of 
22.5 K. The curves are the same as those in Figured] 
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Fig. 3. — Plots of the 1^(240 /zm)/I( 13 CO) ratio versus the atomic hydrogen fraction, 
N(HI)/[N(HI)+2N(H 2 )] are shown for a subsample of the high signal-to-noise positions 
with a 140/zm/240/im color temperature, T dc , near 18 K. Specifically, this sample of points 
is higher than 5a in 1^(140 /im), 1^(240 /im), and I( 13 CO), higher than 3a in I(HI 21cm), 
and with Td c in the range 17 to 19.5 K. The upper plots use the N(H 2 ) values of the non- 
LTE, one-component models and the lower plots use those of the non-LTE, two-component, 
two-subsample models. The left plots include the error bars and right plots omit the error 
bars for clarity. 
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Fig. 4. — These are plots of simulated 1300 fim continuum to 13 CO J = 1 — > spectral line 
ratio data, i.e. 1^(1300 /im)/I( 13 CO), versus simulated 140 /im/240 /xm color temperature 
data, i.e. T^c- The simulations use the parameters listed in the second column of Table 1 
of Paper II and a mass absorption coefficient appropriate for A = 1300 /im (see details 
in Section 14.61) . The upper panel is the plot of 1^(1300 /im)/I( 13 CO) versus Td c for the 
temperature range T dc = 14 to 30 K. The lower panel plots the same quantities, but for the 
larger temperature range of Td c = to 200 K. The lower panel also shows the model curve 
out to T dc = 200 K. 
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Table 1. Best-Fit Parameter Values for the Different Models 



One-Component Models 







LTE 










Full Sample 








AT 




V 






-4 


16.5 


673 








T H > 20 K 








9 


9.0 


140 








LVG 










Full Sample 






AT 


N( 13 CO) 
Av 


n(H 2 ) 




V 


-f 


3 x fO 15 


1 x 10 5 


16.9 


671 






T H > 20 K 






-3 


3 x fO 15 


6 x 10 3 


10.0 


139 



Two-Component Models 



Full Sample 

AT r T N co( 13 CO) „ ' N cl (i3 C Q) 2 

A1 C IdO Ave Uc0 A^c Ucl X " V 

0.04 18 5 x 10 16 2 x 10 4 8 x 10 15 1 x 10 3 5.7 667 

Two Subsample 
Th < 20 K 

AT r T N co( 13 CO) „ N cl ( 13 CO) 2 

A1 C J-dO Ave Uc0 A^c Ucl X " V 

1.0 18 5 x 10 15 1 x 10 5 2 x 10 16 1 x 10 5 4.6 525 

T H > 20 K 

0.4 18 5 x 10 14 1 x 10 4 5 x 10 15 6 x 10 3 8.2 135 



Note. - AT and T d0 are in units of Kelvins. The N( 13 CO)/Av 
quantities are in units of 13 CO molecules ■ cm~ 2 • (km • s -1 ) -1 . The 
n quantities are in units of H 2 molecules • cm -3 . All two-component 
models used the LVG code. See Paper I for discussion of the formal 
and systematic uncertainties. 
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Table 2. Best Estimates of Parameter Value Ranges 1 



Parameter 




Range of Values 


AT b 




-1 to +2K 


Tdo 




16 to 19 K c 


N c0 ( 13 CO) 


2.0 x 10 14 


to 5.0 x 10 15 13 CO cm- 2 ■ (km • s^ 1 )" 1 






> 20 H 2 cm" 3 


N c i( 13 CO) d 


3 x 10 15 


to 2 x 10 16 13 CO cm" 2 • (km • s" 1 )" 1 


n cl 




> few x 10 3 H 2 cm -3 



a See Paper II for details. 

b Assuming two-component models applied to both subsamples. 
c See Section [3TB1 

d For the two-component models applied to the two subsamples, 
the Ncl ^ v CO ' > value would be at the higher end of this range for the 
Tdc < 20 K subsample and at the lower end for the Td c > 20 K 
subsample. 



